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ABSTRACT 

Molecular outflows driven by protostellar cluster members likely impact their surroundings and 
contribute to turbulence, affecting subsequent star formation. The very young Serpens South cluster 
consists of a particularly high density and fraction of protostars, yielding a relevant case study for 
protostellar outflows and their impact on the cluster environment. We combined CO J = 1 — 0 observa¬ 
tions of this region using the Combined Array for Research in Millimeter-wave Astronomy (CARMA) 
and the Institut de Radioastronomie Millimetrique (IRAM) 30 m single dish telescope. The combined 
map allows us to probe CO outflows within the central, most active region at size scales of 0.01 pc 
to 0.8 pc. We account for effects of line opacity and excitation temperature variations by incorpo¬ 
rating 12 CO and 13 CO data for the J = 1 — 0 and J = 3 — 2 transitions (using Atacama Pathfinder 
Experiment and Caltech Submillimeter Observatory observations for the higher CO transitions), and 
we calculate mass, momentum, and energy of the molecular outflows in this region. The outflow 
mass loss rate, force, and luminosity, compared with diagnostics of turbulence and gravity, suggest 
that outflows drive a sufficient amount of energy to sustain turbulence, but not enough energy to 
substantially counter the gravitational potential energy and disrupt the clump. Further, we compare 
Serpens South with the slightly more evolved cluster NGC 1333, and we propose an empirical scenario 
for outflow-cluster interaction at different evolutionary stages. 

Subject headings: ISM: individual objects (Serpens South), jets and outflows, molecules — stars: 
formation, protostars — techniques: interferometric 


1. INTRODUCTION 

1.1. Molecular outflows and protostellar clusters 

Outflows and jets occur during the early stages of pro¬ 
tostellar evolution, and they facilitate the accretion of 
material onto the p rotostar by carrying away excess an¬ 
gular momentum (Bachiller|J1996). A molecular outflow 
consists of ambient gas that is entrained by a protostellar 
jet, expected to occur during the Class 0 and I stages of 
star formation. Class II sources also power jets and winds 
and may drive outflows, however molecular outflows are 
less prominent as a source evolves and the surrounding 
molecular gas becomes depleted. Observationally, out¬ 
flows serve as signposts that emanate from their driv¬ 
ing protostellar sources and indicate that star formation 
is ongoing in a particular region. Further, the charac¬ 
teristics and morphology of outflows may indicate the 
evolution of the forming pr otostar an d constrain the gas 
entrainment mechanism ( Arce et al. 2007). 

Carbon monoxide is the most common probe of molec- 
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ular outflows, because after H 2 it is the most preva¬ 
lent molecule in molecular clouds and its low-level ro- 
tatio nal states are easily excited in these environments 


(e.g., Wilson et al. 2009). The CO molecular outflow 
components are generally characterized as standard high- 
velocity (SHV) with velocities of a few to about 20 
km s _1 , and extremely high velocity (EHV) with even 
hi gher velociti es, typically up to ^ 100 km s -1 (see 


Bachiller 1996). Typically the SHV outflows are more 
massive than EHV, and in this paper we conduct a cen¬ 
sus of mostly SHV components within the protostellar 
cluster Serpens South in order to determine the impact 
of outflows on the surrounding cloud environment. 

Outflows may reach distances up to several parsecs, 
affecting the surrounding protostellar environment cor¬ 
responding to these scales, including the core and cloud. 
Since outflows sweep up ambient gas, they influence the 
intimate relation between a protostar and its surround¬ 
ings. Considering that the majority of clustered star for¬ 
mation occurs within clusters of 100 or more members, 
and the size scales of clusters are typically less than a few 
parsecs (Lada & Lada 2003), there is ample opportunity 
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for outflows to encounter other nearby outflows and their 
respective protostars, as well as potentially disrupt the 
surrounding cloud environment. 

It remains a challenge to quantify the cumulative im¬ 
pact of protostellar outflows on molecular clouds and 
subsequent star formation, but evidence of outflow-cloud 
interaction h as been observed i n a variety of star-forming 
regions (e.g., Arce et al. 2010 Nakamura et al.||2011a|b 


et al. (2008) have assumed a distance to Serpens South 
of 260 ± 37 pc, which is the d istance determined u sing 
photometry of Serpens Main by|Straizys et al.|([l996) and 
consistent with the distance determined for the Serpens- 


Aquila Rift b y Straizys et al. (2003). Most recently 


Dzib 


et al. (20111) used Very Long Baseline Array trigono 


Duarte-Cabral et al. 2012). Specifi cally, outflows interact 


with the ambien t cluster gas ( Fuller fc Ladd||2002 Arce 


& Sargent 2006), i nject momentum and energy into th e 


surrounding cloud ( |Arce et ah 2007 Swift & Welch 2008), 


and they may feed turbule nt motions that aff ect ongo- Dzib et al. (2011). 


metric parallax measurements to determine a distance 
of 429 =b 2 pc to the YSO EC 95 associated with Serpens 
Main, updated from the previous measurement of 415 pc 
by Dzib et al. (2010). Throughout this paper, we assume 
that Serpens South is part of the same complex as Ser¬ 
pens Main , and we adopt the distance of 429 zb 2 pc from 


ing star formation (see also Frank et al. 2014). Models 
of clustered star formation must incorporate outflows in 
order to drive turbulence and produce outcomes close to 
what is observed in star-forming regions, includin g star 


Outfl ows in this region were observed by Nakamura 


et al. (2011a) at an angular resolution of 22', and pre- 


sented in the context of the filamentary 1.1 mm dus t con- 
tinuum image by ASTE. Near-IR observations by |Teix-| 


eira et al. 


mura 2006 |Nakamura & Li 

to 

o 

o 

u2 

ICarroll et al. 2009; 

emission-line objects. Polarization observations by bug- 

Wang et al. 2010: Nakamura < 

WLi 2011 Krumholz et al. 

itani et al. (2011) reveal the importance of the magnetic 


ticularly in clustered star-forming regions, can be con¬ 
strained by observations of clusters. 

In the work presented here, we study outflow charac¬ 
teristics and assess their impact within Serpens South, 
which is the second case study in our survey of clustered 
star-forming regions rangi ng in evolutiona ry stage, be¬ 


ginning with NGC 1333 (Plunkett et al. 2013). NGC 
1333 is a prototypical region that has been well-s tudied, 


while Serpens So uth is more recently discovered (Guter- 
muth et al. 2008), more massive and with a higher frac¬ 


tion of protostars in their earliest stages compared with 
NGC 1333. We describe Serpens South in this context 
in the following sub-section. The goals of the current 
paper are: (a) measure the outflow activity in Serpens 
South, (b) investigate the sustained impact of outflows 
compared with turbulence and gravity, and (c) establish 
an empirical scenario of cluster evolution based on the 
comparison between Serpens South and NGC 1333. 

1.2. Description of the Region: Serpens South 

The protostellar cluster Serpens South merits study for 
its youth, active star formation and corresponding out¬ 
flows, the concentration of young stellar objects (YSOs) 
along a dense filamentary structure, and its relative prox¬ 
imity (429 pc, as discussed in t he following paragraph) . 
Serpens South was discovered by Gutermuth et al. (2008) 
in imaging of the Serpens-Aquila Rift with Spitzer 1RAC, 
in which its dark filamentary structure appeared as ab¬ 
sorption in the 8 /im band. Subsequent studies have 
confirmed the filament ary str ucture of th i s clo u d using 
dust continuum maps (| Andre et al.||20 10l |Maury et al 
201 1 | ) and molecular line emission maps (|Kirk et al.|2013 ' 


Tanak a et al.||2013 Fernandez-Lopez et al.|2014 ). 


Serpens South is one of only a few active star-forming 
regions known within a distance < 1 kpc away. However, 
various distances have been adopted throughout the lit¬ 
erature. Serpens South and the Serpens Main cluster (3° 
to the north) were found to have similar local standard 
of rest velocities ( vlsr ), in addition to their close prox¬ 
imity on the plane of the sky, and therefore the distance 
to Serpens Sou th is based on the distan ce determined for 
Serpens Main (|Gutermuth et al. 2008). The majority of 
works since the discovery of Serpens South by [Gutermuth 


field is well-ordered and oriented perpendicular to this 
filament. 

The classificati on of candidate YSOs in Serpens South 


was first done by Gutermuth et al. (2008) follow i ng the 
data reduction pipeline presented by 
which was developed for the Spitzer 


Evans et al. (2007) 
c2d and Gould Belt 


projects, incorporating 2MASS and Spitzer 1.25 — 70 /im 
phot ometry. Within the 14' x 10' region that they stud¬ 
ied, Gutermuth et al. (2008) reported 54 Class I proto¬ 
stars (including flat-spectrum sources) and 37 Class II 
YSOs. Hence, about 60% of the YSOs are protostars, 
but at the cluster center this fraction rises to about 80%, 


which is where we focus our study (see also Section 4.4). 
Gutermuth et al. ( 2008| ) reported a mean YSO surface 
density of 430 pc -2 at the distance of 260 pc, which cor¬ 
responds to ^ 160 pc -2 at the farther distance of 429 pc 
that we assume here. The median projected distance be- 
tween nearest-neighbor YSOs in this region is 1372 (see 
Gutermuth et al. 2008] for discussion), or 0.03 pc at a 


distance of 429 pc. This nearest-neighbor spacing is less 
than the typical extent of outflows, and therefore this re¬ 
gion provides an ideal case study where outflows likely in¬ 
teract with one another, the surrounding protostars and 
their cluster environment. 


The Herschel Gou ld Belt Survey (Andre et al. 2010 


Bontemps et al. 2010) revealed several of the youngest 


Class 0 members of Serpens S outh, further suggesting 
the overall youth of the region. Maury et al. (2011) sug¬ 
gest that the high number of protostars compared with 
the total number of YSOs in Serpens South indicates 
that this cluster experienced a recent burst of star for¬ 
mation, particularly along a filamentary structure span¬ 
ning from the southeast to northwest. They used com¬ 
plementary 1.2 mm continuum maps obtained with the 
MAx-Planck Millimeter BOlometer array (MAMBO) on 
the Institut de Radioastronomie Millimetrique (IRAM) 
30 m telescope, along with Herschel and Spitzer maps, to 
perform a source extraction and classification over wave¬ 
lengths from 2 /im to 1.2 mm. They determined that this 
cluster has formed 9 Class 0, 48 Class I, and 37 Class II 
objects within an area of 15' x 25' that has a mass of 
~ 610 M©. They suggest that a high SFR and a rel- 
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atively low star formation efficiency are evidence of the 
early evolutionary stage of this region. 

Here we present the highest resolution observations 
published to date of millimeter-wavelength molecular line 
and continuum emission observations made toward the 
most active central region of Serpens South. Specifi¬ 
cally, we observed several isotopologues and energy level 
transitions of carbon monoxide, tracing the SHV outflow 
components. 12 CO traces cool outflow emission, a probe 
of swept up molecular gas. 13 CO is less optically thick, 
and was used to obtain reliable estimates of outflow mass. 
Together, observations of 12 CO and 13 CO provide a com¬ 
prehensive view of the molecular gas kinematics in the re¬ 
gion. Continuum emission was simultaneously observed 
to detect dense dust envelopes around protostars. Some 
of these protostars, along with the protostars observed 
in the previously mentioned studies, may be responsible 
for driving the outflows detected in line emission. 

2. OBSERVATIONS AND DATA 

The methods and results we present here pertain to 
data obtained with both interferometry and single dish 
telescopes, including detections of two energy transitions 
of each molecular isotopologue. In the following sec¬ 
tions we describe the specific observations with the Com¬ 
bined Array for Research in Millimeter-wave Astronomy 
(CARMA), IRAM 30 m, Atacama Pathfinder Experi¬ 
ment (APEX), and Caltech Submillimeter Observatory 
(CSO), and the process of combining interferometer and 
single dish data of the same isotopologue. To put the 
observations into context, in Figure pi w e show the map 
coverages that we describe here; Tablejl] gives a summary 
of the observations. 

We also observed the J = 3 — 2 energy level transition 
of 13 CO with the CSO, and 12 CO with the APEX. Ob¬ 
servations of this higher energy transition are utilized to 
determine velocity-dependent opacity described in Ap¬ 
pendix [A| and excitation temperature according to the 
method described in Appendix |B] 

2.1. CARMA Observations 

Interferometric observations were made with CARMA 
in the D-array configuration during the period 2011 
April and May. We mosaicked a region of approxi¬ 
mately 5' x 1' (oriented with position angle —30°, see 
Figure [l]) within Serpens South. Our mosaic was cen¬ 
tered at R.A.= 18 /l 30 m 00 s , decl.= —2°2 / 0 // and consisted 
of 66 pointings in a hexagonal-packed pattern; the map 
achieved better-than Nyquist sampling with horizontal 
spacing of 30" and vertical spacing of 25"8. In total, we 
observed 33.8 hours, with 24.6 hours on source. Integra¬ 
tions were 30 seconds per mosaic position, and 3 min¬ 
utes on the phase calibrator 1743-038, with 33 pointings 
observed in between each phase calibrator observation. 
Uranus was the primary flux calibrator. 

We simultaneously observed the J — 1 — 0 transitions 
of 12 CO (115.27 GHz) and 13 CO (110.20 GHz) in the 
upper side band (USB). The spectral windows of 12 CO 
and 13 CO had widths of 31 MHz (^ 80 km s -1 ) and 8 
MHz (^ 20 km s -1 ), respectively. With 383 channels 
in each CO window, the velocity resolutions were 0.2 
and 0.06 km s _1 , respectively. In addition, we simul¬ 
taneously observed the continuum emission using eight 



RA (J2000) 


Fig. 1.— C overage of observations, overlaid on map of Spitzer 
8/xm emission (|Gutermuth et al.|2008|). Blue (solid) inner contour 
marks the CARMA mosaic map where Nyquist sampling and best 
sensitivity were achieved, and blue (solid) outer contour marks de¬ 
graded sensitivity at the edge of the map. The inner contour of this 
map comprises the majority of molecular outflow emission. Green 
(dashed) contour marks the IRAM 30m map coverage, which is 
composed of nine 2' x 2' and one 2' x V OTF maps. Yellow 
(dash-dot) contour marks the APEX 12 COJ = 3 — 2 map. Or¬ 
ange (dashed) circles mark four CSO pointings of 13 COJ = 3 — 2 
observations, and diameters correspond to the CSO beam. 


spectral windows of 495 MHz each, resulting in a map 
with a total continuum bandwidth of 4 GHz that covers 
the same region as that of CO. The beamsize of ~ 5" 
corresponds to a physical resolution of 0.01 pc, or 2000 
AU (at a distance of 429 pc). __ 

Data were reduced using MIRIAD (Sault et al. 1995]). 
Anomalous amplitudes and phases were Ragged, as were 
any observations with system temperatures higher than 
500 K. No shadowed data were used, nor were observa¬ 
tions when the source was at elevations between 85° and 
90°. Flux and gain calibrations were applied, and the 
data were Fourier transformed. 

Images of the data were obtained using the task 
invert, specifying a cell size of 2", and then deconvolved 
with the maximum entropy deconvolution mosmem. The 
clean images were then restored with a Gaussian beam. 
We subtracted the continuum image from the molecu¬ 
lar line maps using avmaths in the image domain. A 
summary of the CARMA observations is given in Table 

m 


2.2. IRAM Observations 

We mapped the same region with the IRAM 30 m tele¬ 
scope in 2012 April. The center of the map is R.A.= 
18 /l 30 m 03 s , decl.= -02°02' to match the CARMA map 
coverage. We mapped the region using individual 2' x 2' 
on-the-fly (OTF) maps, with each OTF map scanned 
twice (in orthogonal directions). The 2' x 2' maps cover a 
contiguous region of 38 arcmin 2 outlined in Figure [I] We 
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TABLE 1 

Observations Summary 


Line 


Telescope 

Date 

Rest Frequency 
(GHz) 

HPBW 

maj(") min (") 

Beam PA 
(°) 

Channel Width 
(km s -1 ) 

Bandwidth 

(MHz) 

rms Noise 

12 CO J = 1 

- 0 

CARMA 

2011 Apr, May 

115.27 

5.3 

4.1 

7.4 

0.2 

31 

0.2 Jy beam -1 

13 CO J = 1 

- 0 

CARMA 

2011 Apr, May 

110.20 

5.3 

4.5 

-19.5 

0.06 

8 

0.3 Jy beam -1 

12 CO J = 1 

- 0 

IRAM 

2012 Apr 

115.3 

22.5 

22.5 


0.53 

4000 a 

0.2 K 

13 CO J = 1 

- 0 

IRAM 

2012 Apr 

110.2 

23.5 

23.5 


0.53 

4000 a 

0.05 K 

12 CO J = 3 

- 2 

APEX 

2009 Oct 

345.8 

19.2 

19.2 


0.2 

180 

0.2 K 

13 CO J = 3 

- 2 

CSO 

2011 Jun 

330.6 

35 

35 


0.06 

500 

0.4 K 

2.7mm continuum 

CARMA 

2011 Apr, May 

110.7 

5.3 

4.5 

0.54 


3958 

0.001 Jy beam -1 


a Total bandwidth of the observations. Only 1024 channels centered at = 10 km s ^ were inspected. 


tested several off positions, with the best being located 
at R.A = 18 /l 27 m 3!ll, decl.= -02°22'00" (an offset of 
—2700" in R.A. and +1200" in deck from the central 
map position). 

Here we present d ata obtained using the Eight Mixer 
Receiver (EMIR, |Carter et al. 2012) with the E090 band 
and the backend ETb. We tuned to 110.5 GHz in the 
upper-inner sub-band to detect 13 COJ=l — 0 , and we 
simultaneously observed 12 CO J = 1 - Oin the upper- 
outer sub-band. The velocity resolution was Av = 0.53 
km s -1 , and bandwidth was 4000 MHz (or approxi¬ 
mately 10,000 km s -1 ). The beamsizes at the frequencies 
of 12 CO J = 1 — 0 and 13 CO J = 1 — 0 observed with 
the 30m telescope are 22'/5 and 23 / /5, respectively, corre¬ 
sponding to ~ 0.05 pc at d = 429 pc. 

We reduced the data using the GILDAS/CLASS soft¬ 
ware package (the Grenoble Image and Line Data Anal¬ 
ysis System) 0 For 12 CO J = 1 — 0 we corrected for beam 
efficiency of 77 = 0.78 and forward efficiency of 0.95, and 
for 13 CO J = 1 — 0 we corrected for bea m e ffic iency of 


77 = 0.79 and forward efficiency of 0.95 (Kramer et al 


2013). We fit a first order baseline to all spectra, and 
combined all OTF pointings within 12". The resulting 
rms noise level of the final reduced maps are 0.2 and 
0.05 K per velocity channel for 12 CO and 13 CO maps, 
respectively. 

Deep observations (rms~0.1 K for 12 CO, and 0.04 for 
13 CO) at the off position, using the frequency-switch 
mode, showed that there was in fact significant (^ 2 K, 
with a width of ^ 1 km s -1 for 12 CO) extended emis¬ 
sion in the off position at cloud velocities. We added the 
frequency-switched spectrum from the off-source position 
to correct the OTF maps at these velocities. Although 
the emission at the off position may marginally affect 
the spectra in our map at cloud velocities, at velocities 
beyond ±3 km s -1 with respect to the ambient cloud ve¬ 
locity (the velocities we consider for outflows) we do not 
detect significant emission at the off position. 


2 . 3 . APEX Observations 

The J — 3 — 2 transition of 12 CO was observed with 
APEX in 2009 October. The ^ 6 ' x 6 ' region centered 
at R.A.= 18 /l 30 m 03 s , decl.= — 02 ° 02 ' (see Figures [l] and 
[3| was mapped using the OTF method. The velocity 
resolution of the observation was Av = 0.2 km s -1 , with 
a bandwidth of 180 MHz, or 160 km s -1 . The beamsize 

3 http://www.iram.fr/IRAMFR/GILDAS 


at 345.8 GHz was 19'/2, and the data cube has a pixel 
size of 9". 

For the purpose of obtaining the opacity for each ve¬ 
locity channel (see Appendix [A]), the cube was smoothed 
to have a beamsize of 34'/ 8 , ana regridded to a pixel size 
of 11 " to match the CSO observations (see below). The 
rms of the smoothed data cube is 0.16 K per 0.2 km s -1 
channel. 

These data were also used to determine excitation 
temperatures (see Appendix [B| in conjunction with the 
IRAM 12 CO J = 1 — 0 data, and for this purpose the 
APEX data were smoothed to a beamsize of 2275 and 
regridded to a pixel size of 11 " and velocity resolution 
of Av = 0.53 km s -1 . The resultant rms following this 
smooth and regrid was 0.13 K per 0.53 km s -1 channel. 

2 . 4 . CSO Observations 

The J = 3 — 2 transition of 13 CO was observed 
with CSO in 2011 June. Four positions were ob¬ 
served, with the central pointing being RA=18 /l 30 m 03 s , 
Dec=— 02 ° 02 ' (see Figure HI. The additional three point¬ 
ings had offsets of (+3",-4T'), (-7075,+120"), (-33",- 20 ") 
with respect to the central pointing, selected based on 
various features seen in 12 CO emission. The beamsize 
at 330.6 GHz is 35". Velocity resolution is Av = 0.06 
km s -1 , and the bandwidth was 500 MHz, or 450 km s _1 . 

For the purpose of obtaining the opacity for each veloc¬ 
ity channel, the cube was regridded to match the velocity 
resolution of the APEX observations, Av = 0.2 km s -1 . 
The resulting rms is 0.3 K. 


2 . 5 . Combining Interferometer and Single Dish Data 

Combining interferometer (CARMA) and single dish 
(IRAM 30 m) 12 CO and 13 CO J = 1 — 0 data recovers 
extended emission (up to the combined map size of about 
0.8 pc) that was filtered out by the interferometry, and it 
is imperative for accurately determining mass, momen¬ 
tum and energy of the outflows, presented in Section [+3| 
T he combination follows the same method as described 
Plunkett et al. ( |2013 ), correspo nding to the nonlin - 


ear joint deconvolution method of |Stanimirovic| (| 2002 ). 


This method is particularly useful for recovering flux in 
mosaics that cover regions with clumpy structure. The 
combination was done with MIRIAD. First we regrid¬ 
ded the interferometer map to the velocity resolution of 
the single dish map, and then regridded the pixels of 
the single dish map using the interferometer map as the 
template. 
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Fig. 2.— (a) Serpens South regio n as observed with S vitzer IRAC 8/xm (Gutermuth et al.|2008), overlaid with sources from the literature 
according to their clas sifications by |Maury et al.| ( |2011| shown with star symbols) and m the Spitzer Gould Belt survey (Evans et al. 2007] 
|Gutermuth et al.|2008| shown with circles). The white dashed box marks the region where the majority of CO outflows are lound, shown m 
FigureThe red dashed box shows the region where we detect continuum sources, shown in more detail in the right-hand panel. The 2.7 
m m con tinuum so urces (see Table J2|) are given by crosses and diamonds, such that crosses indicate previously identified YSO candidates 
(Evans et al.|2007| |Maury et al.|2(jTl|, and orange diamonds mark CARMA-2 and CARMA-4 that previo usly were not identified as YSO 
candidates, (b) 2.Y mm continuum (contours) observed with CARMA, overlaid on the Spitzer 24pm map (|Gutermuth et al.|2008|). Noise 
rms of the continuum map is 1 mJy/beam, and the beam is shown in the lower right corner. Contours begin at a level ol 4cr and increment 
by 8cr. Overlaid are the positions of 2.7 mm continuum sources, as described for the left panel. 


The joint deconvolution was done with the MIRIAD 
task mosmem using the dirty interferometer map and sin¬ 
gle dish (regridded) map as inputs, along with dirty 
beams for each. The beam sizes were 5.3" x 4.1" for 
12 CO and 5.3" x 4.5" for 13 CO. We converted the single 
dish maps to have units of Jy beam -1 , in order to match 
the units of the interferometer maps. Several other pa¬ 
rameters were necessary inputs in MIRIAD so that the 
convolution process converged within a reasonable num¬ 
ber of iterations. One such input, “rmsfac”, specifies the 
ratios between the true noise rms (output by MIRIAD 
using imstat in several emission-free regions) and the 
theoretical rms (determined by mossen) for both the in¬ 
terferometer and single dish maps. These factors are 
necessary in order for mosmem to reduce the residuals in 
the combined map (i.e., the difference between the dirty 
image and the model modified by the point-spread func¬ 
tion) to have the same rms as the theoretical rms multi¬ 
plied by the “rmsfac”. We used “rmsfac” of 1.6 and 1.5 
for the interferometer and single dish maps, respectively. 

Another parameter is the flux calibration factor, or 
the number by which the single dish data are multi¬ 
plied in order to convert to the same scale as the in¬ 
terferometer data. We determined this factor using the 
task immerge, choosing uv-data in the range 12 — 25 m 
(which is the range of “baselines” that the interferom¬ 
eter and single dish data have in common), and spec¬ 
ifying channels where emission structure is comparable 
in the interferometer and singe dish maps. We used a 
factor of 1.1 for the 12 CO deconvolution, and 1.3 for the 


13 CO deconvolution. Using these parameters, the decon¬ 
volution converged within 100 iterations for each chan¬ 
nel. The final “CLEAN” maps were generated with the 
task restor. The resulting mean rms sensitivity of the 
combined CARMA and IRAM maps of 12 CO and 13 CO 
were 0.16 Jy/beam/channel and 0.11 Jy/bearn/channel, 
respectively, with 0.5 km s -1 channel width and beam- 
size 5.3" x 4.6" (corresponding to ^ 0.01 pc at d = 429 
pc). The zeroth-order moment map of the combined 
CARMA+IRAM 12 CO data is shown in Figure [ 3 ] (with¬ 
out including the ambient velocity channels, ana where 
the blue and redshifted wing emission are shown in con¬ 
tours), and channel maps of 12 CO and 13 CO are shown 
in Figures [4][5j respectively. 

3. RESULTS 

3.1. Continuum Sources 

Using our CARMA continuum observations, we de¬ 
tected seven distinct continuum sources with flux greater 
than 4cr, and we labeled them in numerical order with 
increasing R.A. The sources are shown in Figure [2] and 
given in Table [2] These sources were fit as 2D Gaussians 
using the MIRIAD task imf it. Using the integrated flux 
measured with imf it, we calculated masses of the con¬ 
tinuum so ur ces ac cording to the following relation (see 
Schnee et al.pQ lO): 

d 2 S u 

B v (T d )k v ’ 


M = 


(1) 
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TABLE 2 

Continuum sources 


Source 

Position 


Size 

PA 

Peak Intensity 

Total Flux 

Mass 

c2d a 

Other 


ct(J2000) 

4(J2000) 

Maj (" 

] Min (") 

(°) 

(mjy/beam) 

(mJy) 

( M o) 

Class. 

Names 

CARMA-1 

18:30:01.2 

-02:03:43.3 

6.6 

2.9 

38.9 

3.6T0.8 

7.1 

0.07 

YSOc 

4 b , Y1 c 

CARMA-2 

18:30:01.4 

-02:01:58.7 

6.5 

3.3 

-61.6 

3.6T1.3 

7.5 

0.07 



CARMA-3 d 

18:30:02.5 

-02:02:47.5 

8.0 

3.2 

18.6 

17.5T2.4 

39.1 

0.36 


SerpS-MM16 e 

CARMA-3 SW d 

18:30:02.5 

-02:02:48.3 

4.0 

3.0 

-87.4 

19.1T2.4 

29.5 

0.27 



CARMA-3 NE d 

18:30:02.6 

-02:02:41.4 

4.0 

3.1 

80.1 

6.4T2.1 

10.0 

0.09 



CARMA-4 

18:30:02.8 

-02:04:08.8 

f 



6.2T1.7 

7.1 

0.07 s 

Galc_red 


CARMA-5 

18:30:03.4 

-02:02:44.3 

5.6 

3.6 

-14.8 

7.0T1.1 

13.1 

0.12 

YSOc_red 

1 b , P2 c 

CARMA-6 

18:30:03.6 

-02:03:08.5 

7.0 

4.6 

36.7 

19.4T1.5 

46.5 

0.43 

red 

SerpS-MM18 h , 

2 h , P3 c 

CARMA-7 

18:30:04.1 

-02:03:02.5 

6.0 

4.1 

-32.9 

63.6T4.4 

131.4 

1.21 

red 

SerpS-MM18 h , 

2 h 

Note. — Position, size, 

position angle 

(PA), peak intensity, 

and total 

flux from Miriad 

task “j 

imfit” with Gaussian 

fit. Size and PA 

are the deconvolved fits. 



[Evans et al.||2007 


From Figure 11 of 

Nakamura et al.|[2011a[, 


pipeline and classification, and Spitzer data are presented by [Gutermuti^et^aJ 


iTeixeira et al.1120121. 


CARMA-3 was fit with a single Gaussian, while CARMA-3 SW and NE correspond to the same systen 
[Maury et al-H201l} classify MM16 as Class 0/1, integrated over a multiple system. 

Deconvolution with Miriad task “imfit” fails. 


fit with a double Gaussian. 


cr 

& We present a mass for this object as if it were an embedded source within Serpens South, although this source was classified a s a, galaxy ca.rj 
d CARMA-6 and CARMA-7 correspond to two MIPS sources, and are identified as a multiple system known as SerpS-MM18 by|Maury et al.| 


Hid ate according to 
|201ll and source 2 


[Evans et al.j200'A 

by|Nakamura et al.| 


where d is distance to the source, S u is the 2.7 mm con¬ 
tinuum flux, and B u (Td) is the Planck function. We 
assumed a dust temperature of Tp = 30 K ( |Chen et al.| 
2013), and a dust opacity of k u = 0.1(V/1200 GHz)^ 
cm 2 g -1 , with /? = 1, whic h corresponds to = 0.0092 
cm 2 g -1 at A = 2.7 mm (Loone y et al.|| 200 QJ). Using 


the dust op acities extrapolated from fOssenko pf fc Hen- 


mng 


( 1994[ Column 6 of Table 1), assuming a gas-to-dust 
ratio of 100 and (3 = 1, then k° h (/3 = 1) = 0.0043 
cm 2 g _1 ; extrapolating instead using /3 = 2.0, then 
k° h (/3 = 2.0) = 0.0021 cm 2 g -1 . In these cases, con¬ 
tinuum masses would be about 2-4 times less than what 
we report, respectively. We note these specifically to re¬ 
mind that without further detailed modeling, our simple 
mass calculation is uncertain by a factor of a few de¬ 
pending on choices of hi and its spectral index /3. |Schnee| 
et al. (2014) found a median emissivity spectral index of 
[3 = 0.9 ± 0.3, and so our assumption of (3 = 1 is within 
their uncertainty. 

Given that our continuum map has an rms of 1 mJy, 
we are able to detect envelopes that have masses greater 
than 0.04 M 0 (at better than a 4cr level). Of the con¬ 
tinuum sources in our map, the sources CARMA-3 and 
CARMA-6/7 were previously detected in_theJVIAMBO 
1.2 mm dust continuum map of Maury et al. (2011), and 
they also appear to coincide with the central clump de¬ 
tected in the ASTE/AzTEC 1.1 mm continuum emis¬ 
sion map (see e.g 


[Kirk et ah]|2013| Tanaka et al.|[2013| ) 

| classified these sources (which they 


Maury et al. (2011’ 
call SerpS-lViM16 and SerpS-MM18) as Class 0/1 and 
Class 0, respectively, based on the envelope mass ver¬ 
sus bolometric luminosity ( M env — Lb n i) evol u tionary di- 


agram that they present. Further, Maury et al. (2011) 


noted that these particular sources are each multiples, 
although they report quantities that are integrated over 
each system. 

The extended morphologies shown in our 2.7 mm con¬ 
tinuum map (see Figure® reinforce the notion that the 


sources CARMA-3 and CARMA-6/7 are each multiples. 
The source CARMA-3 is seen as an elongated struc¬ 
ture with a peak of emission in the southwest and ex¬ 
tended emission to the northeast. According to a double- 
Gaussian fit, the centers of two components in this sys¬ 
tem are separated by ^ 7", or ~ 3000 AU not taking into 
account inclination effects. The total integrated flux of 
the sources is 39 mJy, with the southwest (CARMA-3 
SW) and northeast (CARMA-3 NE) components repre¬ 
senting 75 and 25% of the flux, respectively. These com¬ 
ponents likely share a circumbinary envelope, and since 
there is no absolute minimum in emission marking a clear 
delineation between the two components, in Table [2] we 
also present a single-Gaussian fit and the corresponding 
total mass for this system. A better understanding of 
the individual components merits further modeling. 

For the source CARMA-6/7, although the multiple 
components show some overlapping extended emission, 
we see two distinct peaks of emission, with the north¬ 
eastern peak (CARMA-7) being the strongest continuum 
emission in our map. The peaks of emission are sepa¬ 
rated by ^ 16" (^ 6500 AU) on the plane of the sky, 
but dust emission between the sources may pertain to 
a shared circum-binary envelope. Since we detect two 
distinct peaks, we fit two Gaussians with imfit, and we 
determine the masses of the components CARMA-6 and 
CARMA-7 to be 0.43 and 1.21 M 0 , respectively. 

Two distinct Spitzer MIPS 24 /i m sources are coinci¬ 
dent with the two peaks of mm continuum emission as¬ 
sociated with CARMA-6/7, and the MIP S sources were 
classified as “red” following |Evans et al.| (|2007). The la¬ 
bel “red” indicates that the flux density at 24 /im is at 
least three times the flux de nsity of the neares t avail - 
able IRAC band detection in Gutermuth et al. (2008). 
We suggest that since the MIPS 24 fim sources coincide 
with the continuum sources, they are likely deeply em¬ 
bedded YSO candidates, even though these sources lack 
detection in one of the IRAC bands and therefore are not 
classified as YSO candidates according to the criteria of 
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Evans et al. (2007). CARMA-3 lies about 10" north of 
several Spitzer sources, and at least one of those is par¬ 
ticularly bright in 4 fim and 8 fi m emission, as can be 
seen in Figure |2j Upon close inspection of the location of 
this continuum source, we find weak 4 — 24 /im emission 
at the position of CARMA-3, but it is not clear whether 
the 4 — 24 /im emission is associated with a stronger 
source to the south. This may have detracted from de¬ 
tecting a source at the locatio n of C ARMA-3 according 
to the criteria of Evans et al. (2007), as emission from 
a stronger source could have dominated an underlying 
weaker source nearby. Additionally, a source seen in mm 
continuum emission but not 24 /im emission may signify 
that this source is still IR-dark and therefore at a very 
early stage of evolution. 

Three continuum sources in our map have not been 
previously published as mm continuum sources, but do 
coincide with MIPS 24 /i m sources. These sources 
we name CARMA-1, CARMA-4 and CARMA-5 (see 
Table f2|. CARMA-1 an d CARMA-5 were classified 
as YSOc and YSOc_red (|Evans et all] [2007]), respec¬ 
tively, indicating YSO candidates (“red” is described 
above). CARMA-5 is located only ~ 15" (0.03 pc) from 
CARMA-3, making it difficult to resolve with previous 
single dish observations. Among the continuum sources 
in this region, CARMA-1 has the largest corresponding 
MIPS 24 /im flux. The combination of large 24 /im flux 
and small mass (based on 2.7 mm continuum flux) sug¬ 
gest that this source is among the more evolved mm- 
sources in the map. 

Several of the continuum sources that we have dis¬ 
cussed here have been referred to by various names in 
previous literature, and we refer the reader to the right 
column of Table [2] for these labels. In addition, Teix- 


and CARMA -6/7 - the sources that we detect in com¬ 


eira et al. (2012) indicate an IR source in the region to 
the northeast which they call P4 and they claim drives 
an outflow, but we detect no continuum counterpart for 
this source. To our knowledge, the source CARMA-2 
has not been identified in previous observations. This 
source only has two pixels with continuum flux greater 
than 4 <j, but these are surrounded by pixels with flux of 
2 — 3cr (hence, only marginally detected above the rms 
noise). If this source is indeed a YSO in Serpens South, 
we calculate a mass of 0.07 M©. 

Finally, we also detect unresolved continuum emission 
associated with the source CARMA-4, for which we can¬ 
not confirm whether it is a YSO candidate or a back¬ 
ground galaxy detection. First of all, it was classi fied as 
a galaxy candidate (“Galc_red”) according to Evans et al. 
(2007). This classification was due to its weak MIPS 24 


/im detection of 9.91 4= 1.05 mJy, which is about a fac¬ 
tor of at least 10 weaker than the other MIPS sources 
previously described in this region. However, the coinci¬ 
dent continuum emission may be evidence that this is a 
YSO candidate. If we assume that this source is indeed 
a YSO in Serpens South, the continuum emission implies 
an inner envelope/disk mass of 0.07 M© (see Table |2|. 

In order to compare with sources detected (or not de¬ 
tected) in previous mm- c ontinu um observations of this 
region by Maury et al. (2011|), we scaled the contin¬ 
uum emission fluxes detected at 2.7 mm (110 GHz) to 
the expected flux at 1.2 mm (240 GHz). In that case, 
we expect thermal emission flux to increase by a factor 
of 9 — 20 for /3 = 1 — 2, respectively. For CARMA-3 


mon with Maury et al. (2011) - we measure integrated 
fluxes at 277 mm that, when scaled to 1.2 mm assum¬ 
ing /3 = 1, correspond to 8 0% an d 9 0%, re spectively, of 
the fluxes measured by Maury et al. (2011). This agree¬ 
ment within uncertainties of the 2.7 mm and 1.2 mm 
continuum fluxes, made with interferometer and single 
dish telescope observations, respectively, suggests: (1) 
the majority of emission for these sources arises from 
thermal dust emission on the Rayleigh-Je ans tail (also 
descr ibed in the Orion Molecular Cloud by |Schnee et al. 
2014]), and (2) the sources are compact sources without 
significant extended emission that would be filtered out 
by the interferometer observations. Conversely, assum¬ 
ing /3 = 2 in the scaling from 2.7 mm to 1.2 mm results 
in expected 1.2 mm integrated fluxes that are greater 
by abou t a f actor of two than those reported by |Maury| 


et al. (2011). Without a more detailed study to deter 
mine the actual spectral emissivity of these sources, we 
cannot rule out the possibility that up to half of the 2.7 
mm flux that we detect is contributed from non-thermal 
emission, such as free-free emission. 

We also investigate the remaining sources th at were 
not detected in 1.2 mm observations by |Maury et al.| 
2011). Specifically, we detected the sources CARMA- 
UARMA-2, CARMA-4, and CARMA-5 with a peak 
intensity < 7cr, and we scaled the 2.7 mm fluxes to ex¬ 
pected 1.2 mm fluxes for comparison. Further, account¬ 
ing for beam dilu tion in an 11" beam, the corresponding 
flux detected by Maury et al. (2011) would be reduced by 
at least a factor of about 5 (assuming that an unresolved 
source fully fills our 5" beam, which we consider an upper 
limit for what are likely more compact sources). There¬ 
fore, we find that a la detection in our map corresponds 
to - 2 - 5 mJy/ll"-beam at 1.2 mm, whereas [Maury 


et al. (|201l|) report a la rms noise of 15 mJy/ll"d5eam. 
Therefore7 our 2.7 mm continuum source detections that 


are not in common with those of Maury et al.| (2011) are 
likely due to differences in sensitivity ancfDeamsize of 
the observations; even a 7cr detection in our 2.7 mm map 
would correspond to ^ la detection in a 1.2 mm map 
with an 11" beam. 

In the following section, we refer to several of these 
continuum sources as potential out flow-driving sources, 
and in Figure [3] we show continuum emission alongwith 
0th moment CO maps. In that figure (also Figure |2] and 
M), we overlay the locations of protostars in thejreg ion 
i dent i fied in the Spitzer Gould Belt survey ( |Evans et al. 


2007 Gutermuth et al. 2008). We note that while we 


focus on the continuum sources in the context of our 
discussion of possible outflow-driving sources, there are 
many (^ 45) densely packed protostars in the mapped 
region. Specifically, within a radius of 30" from a point 
centered between the sources CARMA-3, CARMA-5, 
CARMA-6, and CARMA-7 (the location that appears to 
be the approximate apex of the majority of outflow emis¬ 
sion) reside 10 protostars identified with Spitzer , and 18 
protostars reside within a radius of 1'. Whether some 
Spitzer sources correspond to weak continuum sources 
below our detection limit, or do not correspond to con¬ 
tinuum sources at all, they noneth eless may contribute to 
the outflow activity described in §3.21 In a future study, 
we plan to present observations with higher resolution 
and sensitivity in the central region that will likely give 
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Fig. 3.— Integrated intensity maps of the blue and redshifted outflow emission in Serpens South traced by (a) 12 CO J = 1 — 0 emission 
observed with CARMA+IRAM, and (b) 12 CO J = 3 — 2 emission observed with APEX, with beam sizes for each given in lower left corners. 
The scale bar applies for both panels. Grayscale shows integrated intensity with 3 < |u 0 ut| < 31 km s —1 , and blue and red contours show 
the blue and redshifted moment zero maps within this velocity range. Contours begin at 5cr for each respective map, and increase in 
intervals of 10cr. Green contours show 2.7 mm continuum observed with CARMA. beginning with 4a and inc reasing by 8cr. Crosses mark 
the 45 protostars in the region from the Spitzer Gould Belt survey (|Evans et al.|2007 Gutermuth et al.|2008|), and these are also shown in 
Figure [2] according to their classifications. 


us a better grasp of specific out flow-driving sources, and 
in that case we plan to discuss the Spitzer and continuum 
sources further. 

3.2. CO Outflow Morphologies 

Our suite of high spatial-resolution mm-wavelength 
12 CO and continuum maps, along with archival Spitzer 
IRAC images, allow us to scrutinize the outflow mor¬ 
phologies in the most dense regions, assess the impact of 
these outflows on their surrounding environment, suggest 
probable driving sources where possible, and compare 
with previous identifications from the literature. Along 
with the dense concentration of YSOs in the central clus¬ 
ter previously mentioned, the resolution of our observa¬ 
tions and the breadth of the outflow emission make it a 
challenge to identify individual driving sources. In this 
section we describe the CO outflow morphologies that we 
observed at particular velocity ranges, and we present the 
qua ntitative characteristics of these outflows in Section 

We accompany the following discussion with references 
to several CO maps in Figures [3][9j beginning with Figure 
[3] showing 0th moment maps oF 12 CO J — 1 — 0 and 
J = 3 — 2 emission. In Figure [ 4 J we show 12 CO J = 1 — 
0 velocity-binned maps, summing channels in variable¬ 
sized bins so that near the cloud velocity we sum over 
two consecutive channels (Av = 1.1 km s -1 ), and for 
high outflow velocities we sum over ten channels (Av = 
5.8 km s _1 ). Additionally, 13 CO emission is shown in a 
similar manner in Figure [5j while featuring especially the 
lower-velocity channels where 13 CO emission is detected. 
Figure [6] shows 13 CO overlaid with 12 CO contours for 
certain velocity ranges. Figures [Tp] are moment maps of 


low-, mid-, and high-velocity outflow emission probed by 
12 CO J = 1 — 0 CARMA+IRAM observations. Velocity 
ranges shown in those maps are asymmetric about the 
cloud velocity in order to reveal particular features that 
were identified by close inspection of the data cube, and 
those features are discussed specifically in the text in the 
following sections. An outflow candidate is shown in a 
similar way (but at low er outflow velocities) in Figure [9j 
and discussed in §3.2.3| 

Throughout, we use the notation of vlsr to indicate 
channel velocity (relative to local standard of rest), v c i ou d 
for cloud velocity, and v out = vlsr — ^cloud for outflow 
velocity (along the line of sight). We adopt a systemic 
cloud velocity of v c i ou d = 8 km s -1 per taining to the re¬ 


gion we mapped. Previo us works (e.g., Nakamura et al. 


2011a Kirk et al. 2013) adopted a cloud-wide vel ocity 
centroid of v r inud = 7.5 km s -1 , but Figure 4 of 


Kirk 


et al. (2013) shows that in the region we mapped (cen¬ 
tered at R.A.= 18 h 30 m , decl.= —2°2'), there is a pre¬ 
dominance of a larger centroid velocity of approximately 
Vdoud = 8 km s -1 . 


3 . 2 . 1 . Blueshifted outflows 

Prevalent blueshifted CO features appear north¬ 
east and south of the continuum sources CARMA- 
3/CARMA-5 and CARMA-6/CARMA-7, shown in Fig¬ 
ures |3]|4| and the upper row of Figures |7p| In the 
lower blue velocity channels (i.e., especially with approx¬ 
imately v out = —6 to —4 km s -1 ), these features appear 
as wide-angle outflows with opening angles of ^70° and 
45°, and position angles of ~45° and 180°, respectively. 
However, they may in fact be several more collimated 
outflows super-imposed on one another (see Figures [4j 
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Fig. 4. — 12 CO J = 1 — 0 channel maps using combined CARMA and IRAM data, with positions of 2.7 mm continuum sources observed 
with CARMA marked with plus signs. Velocity channels are grouped such that panels in top and bottom rows have velocity ranges of 
Av = 5.8 km s —1 (except for upper left and lower right, which have Av = 5.3 km s —1 ), and panels in middle two rows have velocity ranges 
of Av = 1.1 km s -1 . Color scale bars are given for each corresponding row. Specific outflow velocities are indicated in respective panels, 
and velocity channels within ±3 km s —1 of the cloud velocity are omitted. 


[7] and |8|. These outflowing structures correspond ap¬ 
proximately to those i dentified as Bl, B2, B3 and B4 by 
|Nakamura et al.| (|2011a|) , which we have also labeled in 
Figure [ZhT Tor reference. The outflow features in low- 
velocity channels suggest that Bl and B2 comprise the 
cavity walls of a single northeastern outflow lobe, and 
similarly B3 and B4 comprise the walls of a southern out¬ 
flow lobe. Wide-angle, far-reaching SHV outflows with 
cavity walls that are nearly parallel at larger distances 
from the driving sources have been observed in other re- 


et al. 2013 

,HH 

(|Lee & Ho 

2005 


2005|) . In the same direction as the outflows 


13 CO emission correspond to the strongest 12 CO emis¬ 
sion, shown in Figure p[a) (see also Figure p|, indicating 
entrained gas even at lower velocities than vmat is probed 
by the optically thick 12 CO. We note also that since pro¬ 
tostars in Serpens South are in such close proximity, and 
given that we expect each protostellar source to drive a 
bipolar outflow, it is likely that what appear as broad or 
wide-angle outflows are in fact comprised of several more 
collimated outflows from the distinct continuum sources 
mentioned here or from the other ^ 10 — 20 protostars 
that reside in the central region but remain undetected 
in continuum observations. 

The northern wall of the blue northeast outflow struc¬ 
ture is comprised of clumpy emission (features A and 
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Fig. 5. — 13 CO J = 1 — 0 channel maps using combined CARMA and IRAM data (with plus signs as in Figure^). Velocity channels 
are grouped such that panels in top and bottom rows have velocity ranges of Av = 1.6 km s -1 (except for upper left and lower right), and 
panels in middle two rows have velocity ranges of Av = 1.1 km s — 1 . Color scale bars are given for each corresponding row. Specific outflow 
velocities are indicated in respective panels, and velocity channels within ±0.6 km s -1 of the cloud velocity are omitted. Additionally, 
dashed boxes indicate the regions shown in more detail in Figure [6] 


B in Figure |7jb)-(c)) along an axis with position angle 
of ~ 20°, appearing to originate from near the sources 
CARMA-3 and CARMA-5 and extend ~ 120" (0.25 pc) 
northeast. At blue velocities with \v out \ > 17 km s -1 , the 
outflow is comprised of delineated clumpy features that 
form part of the structure’s wall. Another clumpy fea¬ 
ture (C in Figure mb)) is seen up to \v out \ = 15 km s -1 
at a location 45" (0.1 pc) northeast of CARMA-3 and 
CARMA-5 (along an axis that intersects CARMA-5 with 
position angle ^ 60°). Outflow emission in the north¬ 
east is seen with \v out \ up to ~ 30 km s -1 (Figure [7^-c), 
whereas the southern blue outflow shows emission up to 
\v 0 ut\ ~ 10 km s” 1 (Figure [7^ a)). 

These clumpy features may be evidence of episodic 
ejection from one or more sources in the central region, 


but we cannot be certain of whether these come from 
one or more driving sources in that region. We also note 
that the clumpy and bullet-like blueshifted emission in 
high-velocity channels is sufficiently strong and localized 
that negative side lobes appear at their edges due to the 
synthesized beam, as seen in the upper left panels of Fig¬ 
ure [4] These negative sidelobes have absolute values that 
are less than 4 sigma in individual channels, and they are 
not included in any cal culations of mass, momentum or 
energy (see Section 3.3). 

Several “infrared H- H objects” (e.g., described by 
Naka mura et al.||2011a ) and H 2 features (Teixeira et al 


2(J12| ) may be associated with t he CO outflows discussed 


here. Nakamura et al. (2011a) report the infrared H-H 
objects, likely resulting from shocks by protostellar out- 
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Fig. 6. — 12 CO and 13 CO emission in several regions of the maps, showing where the low-velocity entrained gas traced by 13 CO coincides 
with molecular outflows traced by 12 CO. Panels (a) and (b) show blue and redshifted 12 CO emission (contours) on pixel maps of 13 CO 
at lower velocities of \v ou t\ < 10 km s _1 . Panel (c) shows redshifted 12 CO (contours) overlaid on lower-velocity blueshifted 13 CO. The 
contours are the same as in Figure [3ja), and regions correspond to dashed boxes shown in Figure [ 5 ] 


flows, which they call K1 and K2 toward the southwest 
of the map. These H-H objects have bow shapes indi¬ 
cating that they are moving toward the southwest and 
probably originated near the regi on of high protostellar 
density at the center of our map. Teixeira et al. ([2012) 


detected two H 2 features (MH03249b and 1V1HO3250) in 
which the emission appears to correspond wit h the IR 
H-H o bj ects K1 and K2 , respe ctively (Nakamura et ah] 


2011a). 


Teixeira et al. 


suggest that F4 drives 
MH(J^249b, Y1 (our cAr 1V1A-1) drives MHO3250 (de¬ 
scribed below), and P3 (our CARMA-6) drives H 2 fea¬ 
tures southward (bey ond the edge of our map) including 
MH03248 (see £3.1 for continuum source labels). This 
implies that the western component of the southern blue 
outflow we detect is driven by CARMA-6 (P3) and the 
easter n component may be related to the source P4 from 
Teixeira et al. (2012). However, in this mapped region 
there are several plausible driving sources, and higher- 
resolution data are needed to determine their origin(s), 
as well as whether these components may form cavity 
walls of the same outflow. Figure [7] marks these features 
from the previous literature with magenta dashed lines 
and ellipses. 

One 12 CO emission feature that appears very com¬ 
pact and is apparent thanks to our high-resolution data 
is located approximately equidistant from the contin¬ 
uum sources CARMA-3, CARMA-5, CARMA-6 and 
CARMA-7 (feature Db in Figure [7^), but not coinci¬ 
dent with any single source detected here or identified 
in the literature to our knowledge. Molecular CO emis¬ 
sion is seen here in velocity channels up to \v out \ = 20 
km s _1 . This outflow’s redshifted counterpart (feature 
Dr, discussed in the following section), is coincident (on 
the plane of the sky) with this blueshifted feature, and 
therefore this may be a bipolar outflow oriented nearly 
along the line of sight. The compact morph ology is sim- 
ilar to that of IRAS 4B in NGC 1333 (Yildiz et al.||2012 
Plunkett et al.||20l3 ). 


3 . 2 . 2 . Redshifted outflows 

Most redshifted outflow emission is located north- 
northwest of the continuum sources in our map (see Fig¬ 
ures [3][4] and lower panels of Figures 7][8), comprising a 


structure with position angle of —20° to —30° and reach¬ 
ing about 210" (0.4 pc) from the continuum sources. At 
the farthest extent of this emission, especially at low ve¬ 
locities with \v out \ < 8.5 km s -1 , the emission has bow¬ 
shaped morphology, and at higher velocities the emission 
is more clumpy (labeled E in Figures [7je) and (f)). This 
far-re aching redshifted emissio n appears to correspond to 
R2 in Nakamura et al. (|2011a|), but the features are seen 
in more detail here. Additionally, in Figure 6] we show 
the expansive 12 CO redshifted emission overlaid on the 
low-velocity 13 CO blueshifted emission. It appears that 
the dense gas traced by 13 CO emission at low-velocities 
is being entrained and spanning an even larger region 
than can be seen with 12 CO (due to its higher optical 
thickness), with especially strong 13 CO emission seen to 
the northwest and south. 

At velocities up to about \v out \ = 12 km s -1 (shown in 
Figures|7[d) and (e), coinciding with the label R3), 12 CO 
emissiontorms an oblong cavity that extends from the lo¬ 
cation of the continuum sources northwest about ^ 100" 
(0.21 pc). It is not clear where this SHV emission orig¬ 
inates, or whether it is the superposition of several out¬ 
flows originating from distinct locations. In low-velocity 
channels, the 13 CO emission corresponds to the peaks 
of 12 CO emission that outline this feature (see Figure 

t ). At higher velocities, the 12 CO structure appears 
impier, while at lower velocities the structure is more 
widespread and diffuse. 

At low red velocities (see Figure [7^d)), we detect a 
very collimated feature with position angle ^30° (i.e., 
northeast) that extends about 140" (0.28 pc) from the 
central cluster of continuum sources and reaches a veloc¬ 
ity of \v nu t\ = 8 km s -1 . This feature is labeled as R1 by 


Na kamura et al.| (2011a ) and coincides with MH03249r 
By" Teixeira et al. |T|M§. A line drawn through the fea¬ 
ture "with approximately the same position angle passes 
close to CARMA-6 and CARMA-7 to the southwest, sug¬ 
gesting that one of these sources (or another one of the 
~ 10 nearby protostars seen with Spitzer ) is a possi¬ 
ble driving source. Another, shorter emission feature to 
the west/so uthwest of MHQ 3249r and corresponding to 
MH03247 
from near 


<1 


Teixeira et al. 


2012), appears to originate 


ARMA-5 and CARMA-3 and reaches a dis- 
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Fig. 7. — 12 CO outflow emission, showing low-velocity (left), mid-velocity (middle) and high-velocity (right) outflowing gas features. 
Blue and red contours represent blue and redshifted gas, and velocity intervals are indicated for each panel, where outflow velocity, v ou t, is 
relative to the central cloud velocity of 8 km s —1 . Contours begin at 5a for each 0th moment map, and increment by 15cr. The beamsize 
and scale bar are given in the lower right corner. Background is Spitzer S/am. emission (Gutermuth et al.|2008|. Boxes 1, 2, and 3 refer 
to the regions shown in Figures [8{a-c),[8{d-f), and [9] respectively. Ellipses and the arrow with solid cyan and pink lines indicate outflow 
features mentioned in the text, ana plus signs and diamonds mark the position s of continuum sources detected by our CARMA observations, 
labeled in panel (c). K1 and K2 are IR H-H objects identified by[Nakamura ~ al.I (2011a). and correspond to the southwestern tips of the 
molecular hydrogen objects MH03249b and MHO3250 identified by Teixeira et al. 2012 g P4 is an IR source. MHO features from Teixeira 
|et al.| |2012~ see their Figure 1) are indicated with magenta dashed lines and ellipses. 


tance of 60"(0.12 pc) northeast of CARMA-5 (see ma¬ 
genta annotation in Figures [8jd)-(f)). The position angle 
of this feature is slightly less than that of the longer colli¬ 
mated flow previously mentioned in this paragraph, and 
it is also possible that these two features form parts of 
the same flow, in which case the flow may be precessing. 

A high-intensity clump is seen equidistant from the 
continuum sources and coincident with the blueshifted 
clump previously mentioned, with redshifted outflow 
emission up to \v out \ = 18 km s -1 (marked as Dr in 
Figures [7je) and (f)). Since we detect emission at this 
location m blue and redshifted velocity channels with ap¬ 
proximately v out = ±20 km s _1 , this appears to be an 


outflow oriented along the line of sight (see Figures |8|d)- 
(f). Figure [7jb) marks the blueshifted counterpart Db 
at the corresponding location on the plane of the sky. 
Although we do not confirm a driving source with con¬ 
tinuum emission, this emission is nearby several MIPS 
24 /i m sources, and hence the outflow may be associated 
with a source that is very embedded. 

3.2.3. Outflow emission near CARMA-1 

Given the relative isolation of CARMA-1 in this field 
(compared with the sources in the center of the clus¬ 
ter, including CARMA-3, CARMA-5, CARMA-6, and 
CARMA-7), we made a close channel-by-channel inspec- 
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Fig. 8. — 12 CO outflow emission, showing low-velocity (left), mid-velocity (middle) and high-velocity (right) outflowing gas features 
within Box 1 (panels a-c) and Box 2 (d-f), respectively (see Figure [7]). Contours and background images are the same as in Figure [7] 
White plus signs show locations of continuum sources, while black crosses indicate Spitzer protostars, as in Figure [3l The beamsize ana 
scale bar are given in the lower right corner. Lines and ellips es indicate outflow features mentioned in the text. The dashed magenta 
line marking MHO 3247 replicates that of|Teixeira et al. (|2012| see their Figure 1), and coincide s specifically to the 8 um emission feature 
slightly southeast of the line. A, B, C, Db, andTJr are outflow features, and P4 is an IR source iTeixeira et al. (2012) that we describe in 
Section |3] 


tion of emission in the immediate vicinity of CARMA-1 
to identify outflows uniquely associated with this Class I 
source. This region is shown in more detail in Figure [9j 
and we highlight the low-velocity channels where outflow 
morpholo gy is appare nt . Prev iously, this source (labeled 


as Y1 by Teixeira et al. 


was proposed to drive the 


H 2 feature 1V1HU32 5U toward the southwest, and Naka¬ 


mura et al. (2011a) labeled two outflow features (R4 and 
B5) that were possibly driven in that direction by this 
source. Qualitatively our data support the previously 
proposed source and outflow identifications, but now we 
detect emission closer to the source, and it suggests that 
there may be additional outflow components associated 
with this source. 

Several clumpy 12 CO redshifted features appear to be 
concentrated near CARMA-1, shown in Figures|9](d)-(f), 
with features typically appearing in only severalcontigu- 
ous velocity channels. For example, at redshifted veloc¬ 
ities \v ou t\ < 5 km s _1 , a narrow feature extends south¬ 
west ^ 35" with a position angle of 240°; this feature 
is oriented along the same direction as the feature R4 
of 


Nakamura et al. 


(2011a) previously mentioned. Ad¬ 
ditionally, blueshifted emission is seen super-imposed in 
low-velocity channels (\v out \ < 6 km s -1 , Figure [9^a)), 
which entails that these features trace the cavity walls 
and this outflow is approximately in the plane of the sky. 
The low radial-velocities that we detect here and the pro¬ 
posed orientation in the plane of the sky are consistent 
with the related MHO being located such that the out¬ 


flow has a relatively large spatial extent 90"). How¬ 
ever, the blueshifted emission feature is oriented slightly 
more south than the redshifted feature, and therefore 
they may not necessarily form the same coherent struc¬ 
ture. 

At blueshifted velocities with \v out \ < 10 km s -1 we 
detect clumpy emission with position angle ^ 25° and 
extending ~ 30" (0.06 pc) northeast from CARMA-1, 
marked as feature F in Figure |9|a)-(c). Emission in this 
direction is also seen super-imposed in redshifted veloc¬ 
ities with \v ou t\ > 7 km s -1 (Figure J9F). At slightly 
higher velocities - seen particularly inFigure [9](b)-(c) 
and marked as feature G - blueshifted molecular out¬ 
flow emission is seen to extend ~ 35" north of CARMA- 
1, with position angle of 0° to —10°. It is not clear 
whether the components with position angles of ~ 25° 
and ~ —10° to the north of CARMA-1 correspond to the 
same outflow with opening angle of 35°, or whether these 
are distinct outflows emanating from multiple sources in 
the proximity. Also, in redshifted velocities shown in Fig¬ 
ure p\e) emission is seen to the north of CARMA-1 with 
a similar position angle, perhaps constituting the far-side 
wall of the same outflow and suggesting that the incli¬ 
nation of this outflow is nearly in the plane of the sky. 
In the 13 CO map, we inspected low-velocity channels in 
the vicinity of CARMA-1 and we found 13 CO emission, 
in particular, to the north of CARMA-1 (see Figure [6]), 
although the emission is somewhat diffuse and does not 
conclusively point to this source. 
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Fig. 9. — 12 CO emission near the candidate outflow from CARMA-1 (cross), corresponding to Box 3 in Figure [7] Contours are in units 
of rms for each moment map, beginning at a level of 4a and incrementing by 4a at low-velocity, and 3a for other panels. The beamsize 
and scale bar are given in the lower right corner. Velocity intervals are indicated for each panel, where ou tflow velocity. v ou t , is relative to 
the central cloud velocity of 8 km s -1 . Features marked with letters and arrows are addressed in Section [3. 2. 3| 


In summary, here we provide additional evidence that 
one or more molecular outflows are being driving from 
the location of CARMA-1. However, a clear picture of 
a bipolar outflow is not immediately apparent. Further, 
the collimated low-velocity blueshifted feature F which 
is a candidate outflow lobe does not have an equally col¬ 
limated counterpart in the red. It seems that the outflow 
structure of this source is either irregular, or the emission 
is being overwhelmed by stronger emission from the sur¬ 
rounding region which makes the outflow appear more 
clumpy rather than collimated. Alternatively, if gas has 
already been swept away from the immediate region by 
other stronger outflows such as those to the northeast, 
there remains little gas to be entrained by this driving 
source CARMA-1. 


depth, T 7 


mb, 12’ 


given by Equation 


A4 


Appendix [B] 


de¬ 


scribes the method(s) to solve for excitation temperature, 
T e x- In the above equations, J = 0 is the lower energy 
level of the rotational transition, and gj = 2 J + 1 = 1 is 
the degeneracy of the lower state. For Xco w e as sume 
a stan dard abundance ratio [ 12 CO/H 2 ] = 10 -4 (Frerking 
et al. |1982 ). For 12 CO J = 1 — 0, v = 115.27 GHz, and 


/i is the 12 CO dipole moment with a value 0.11 Debye. 
Q(Tex) is the partition function which here is approxi¬ 
mated as Q(Tex) = &T e x/hF>, and B = 57.635 GHz is 


the rotation constant for CO. Ej +1 = 7.64 x 10“ 16 erg 
is the energy of the upper level such that the energy of 
the transition is E = Tj + i — Ej, and Ej+i/k = 5.53 K. 

Having measured column density for each pixel, mass 
is then determined to be: 


3.3. Mass, momentum and energy of outflowing gas 

The mass, momentum and energy of outflowing gas 
were calculated in order to assess the overall contribu¬ 
tion of outflow activity within the mapped region, and 
we compared these with turbulent momentum and en¬ 
ergy and gravitational potential energy. First, we calcu¬ 
lated H 2 column density ( Nh 2 ) for each pixel (x,y) and 
velocity channel (with channel width Av) as a function of 
excitation tempera ture (Tex) according to the following 
relations (following |Dunham et al.||2014 ): 


M(x, y, v ) = 2.72 m H N H2 (x, y, v)A, 


(4) 


Nh 2 (x, y, v , T ex ) = /(J, T ex , X CO )T^ M2 Au 

^ 3k 2J + lQ(T ex ) 

j(T, T ex , X co ) — X co l2 87r 3^ 2 J+1 gL 

The opacity correction is explained in Appendix [A[ with 
the main-beam temperature of 12 CO corrected for optical 


where mn is the mass of a hydrogen atom, such_that 
2.7 2 mu is the mean molecular weight of the gas (Bourke 
et al.|19 97), and A is the area of a pixel. We sum over the 
cube m spatial and spectral dimensions including pixels 
that have signal greater than 3cr to get the total outflow 
mass, momentum and energy: 


1 E vei M(x,y,v) (5) 

1 E ve iM(x,y,v)\v - vcloud | (6) 

''ea'EvelM (x,y \v V c loud\ i 



M out = 

(2) 

Pout — 

E J+, 1 . 

Eout 

e 3) 


(7) 


where v is the velocity corresponding to a given chan¬ 
nel, and v c i ou d = 8 km s -1 is the cloud velocity. Mass, 
momentum and energy for velocity channels with 3 < 
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TABLE 3 

Outflow mass, momentum, energy 


Method 

All (Mg) 

Mass 

Blue a (Mg) 

Red 0 (Mg) 

Momentum 
(M® km s -1 ) 

Energy 
(xlO 45 erg) 

(1) Tex(%,y,v), pixel-by-pixel + channel-by-channel c 

2.8 

1.2 

1.6 

13.6 

1.0 

(2) T e x (v), channel-by-channel c 

2.9 

1.2 

1.7 

13.9 

1.0 

(3) T ex = 10 K 

2.2 

1.1 

1.2 

10.8 

0.8 

(4) T ex = 30 K 

4.6 

2.3 

2.4 

22.4 

1.7 

(5) T ex = 50 K 

7.2 

3.5 

3.7 

34.7 

2.6 


a Blue channels correspond to \v ou t\ — 3.3 to 30.9 km s 1 


b Red channels correspond to \v ou t\ = 3.0 to 30.6 km s 1 . 
c Means: Tex(x,y,v) = 17 d= 8 K for (1), and Tex(u) = 15 db 4 K for (2). 


\v ou t\ < 31 km s _1 and the region we mapped cover¬ 
ing ^ 45 arcmin 2 (see Figure flj) are given in Table 
[3] The lower-limit of this velocity range was chosen 
by inspecting low-velocity channels by eye to exclude 
those channels dominated by ambient cloud emission, 
and the upper-limit was dictated by the spectral cov¬ 
erage of our observations. We measure a total outflow 
mass of 2.8 M®, of which blueshifted and redshifted chan¬ 
nels contribute 43% and 57%, respectively. These masses 
were calculated using the Tex(#, 2 /, v) values determined 
with the pixel-by-pixel ± channel-by-channel (velocity- 
dependent) method as described in Appendix p] (a.k.a. 
method 1). Using excitation temperatures determined 
according to the alternative methods 2-5 (these methods 
either spatially average over pixels in each channel, or 
assume a constant excitation temperature, as described 
in Appendix [B| result in total outflow masses ranging 
from 2.2 to 7 TM 0 , listed in Table [3] 

In Figure [TO] we show the outflow mass distribution per 
velocity bin for Serpens South (panels a and b). For the 
blue and redshifted mass spectra, we fit power laws of 
the form m(v) oc v ~ 7 for the velocity range 3 < \v out \ < 
15 km s -1 , and we found slopes of 7 = 3.5 ±0.1 and 7 = 
3.3 ±0.1, respectively. The slope may be an evolutionary 
indicator of mass ej ection i n a clo ud, such that the slope 
steepens with time (Smith et al. 1997). We performed a 
similar analysis of the outflow m ass within the slightly 
more evolved region NGC 1333 (Plunkett et al. 2013), 
assuming T e x = 17 K which is the mean 1 ex in Serpens 
South, and the blue and redshifted mass spectra for that 
region can be fit with power laws with slopes of 7 = 
5.1 ± 0.2 and 7 = 4.9 ± 0.1 (see Figure [TO] panel a and c). 
These two regions are therefore presented as a examples 
for which the shallower (steeper) slope pertains to the 
less (more) evolved region. 

It has been previously suggested that a break in 7 at 
~ 10 km s -1 can be inter preted as two distinct outflow 


spectra are apparent at outflow velocities of ^ 5 — 10 


veloc ity components (e.g., Richer et al. 2000 Su et al 
|2004|). In that case, the high-velocity gas corresponds 
to a recently accelerated component and steeper slope 
( 7^3 — 4), and the low-velocity gas is a slower, coasting 
component with a less steep slope ( 7 ^ 2 ). Neither Ser¬ 
pens South nor NGC 1333 showed a break in their power 
spectra according to our primary analysis. We note that 
our analysis included a velocity-dependent correction for 
opacity of the 12 CO line, which especially accounted for 
more “hidden” mass at low-velocities than in the case of 
the optically thin assumption. When assuming that the 
12 CO line is optically thin (r << 1 ), breaks in the power 


km s 1 (see Figure 10 b-c). Our findings are in line with 


the suggestion by Bally et al. (1999) that a broken power- 


law index for outflows may not be a physical distinction 
but rather an artifact of the optically thin assumption. 

Comparing blue and redshifted mass at each outflow 
velocity, the ratio M out:b iue/M out:re d{\vout\) ranges from 
0.1 to 1.5, with a mean of 0.5. While we expect out¬ 
flows to be bipolar, and hence the blue and redshifted 
gas should total approximately equal amounts, in this 
map it appears that slightly more gas is being entrained 
by redshifted lobes than by blueshifted lobes, especially 
in velocity channels with \v ou t\ ^5 — 10 km s -1 . Possibly 
some outflow gas has escaped the cloud more readily in 
the foreground than in the background, or extends be¬ 
yond the map edges, depending on the locations of the 
driving protostars within their nascent gas cloud. This 
mass ratio may also entail that the surrounding gas en¬ 
vironment was not homogeneous in mass and/or density 
at the time when protostars began to drive outflows. We 
present maps of T e x in Appendix [Bj and the highest T e x 
values are found in redshifted velocity channels. 

In addition to the mass calculation, we calculate out¬ 
flow momentum to be P out = 13.6 M® km s -1 , and en¬ 
ergy to be E out ~ 1.0 x 10 45 erg. Table [ 3 ] lists these di¬ 
agnostics using each T e x method, and we consider these 
values to be lower limits because they are uncorrected 
for inclination angle of the outflows. Assuming an av- 


erage inclination angle with respect to the line of sight 

£ = 57.3° (e.g., Nakamura et al. 

2011a Dunham et al. 

2014), and assuming that all outfi 

low motion is along the 


jet axis (see |Downes fc Cabrit 2007), results in greater 


momentum and energy by factors of 1.9 and 3.4, respec 
tively. Hence, momentum and energy from method 1 are 
Pout = 25.3 M® km s -1 and E out = 3.3 x 10 45 erg, after 
correction for an average inclination angle of £ = 57.3°. 

We performed a comparison of the cumulativ e mass 
an d ener gy of outflows with those determined by |Naka- 


mura et al. ( 2011 a) using 12 CO J = 3 — 2 assuming: ( 1 ) 


constant T e x=30 K, ( 2 ) distance of 260 pc, and (3) in¬ 
cluding only the velocity channels between —15 km s -1 
to 4 km s -1 , and 11 to 30 km s -1 , with a cloud velocity 
v c = 7.5 km s -1 . We found the cumulative mass and 
energy of ou tflows to be 1.8 a nd 1.1 times the values 
presented by Nakamura et al. (2011a), respectively. In 
addition to using a different energy level transition, this 
discrepancy may also be attributable to the_co rrectio n 
we mak e for opacity of CO, whereas Nakam ura et al. 
( 2011 a) assume 12 CO J = 3 — 2 is optically"thin. On tKe 
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Fig. 10.— Outflow mass spectra for Serpens South (circles) and NGC 1333 (triangles), with blue and redshifted velocity channels plotted 
independently with their respective colors. Panel (a) shows both regions for comparison, and aqua and magenta lines (solid for Serpens 
South, dashed for NGC 1333) are the best-fit power-laws to the blue and redshifted mass spectra, respectively. Panels (b) and (c) show 
Serpens South and NGC 1333 separately, where mass spectra with no velocity-dependent correction for opacity (i.e., assuming r << 1) are 
also given by plus and cross symbols for Serpens South and NGC 1333, respectively. 


other hand, we also note that the map coverages are not 
identica l, and specifically the bi ne features B7, B8, and 
B9 from Nakamura et al. (2011a) are not covered entirely 
by our map. Nonetheless, we consider these masses and 
energies to agree within the uncertainties that we cannot 
account for here. 



culation potentially omits some outflowing gas. Making 
an admittedly simple extrapolation to velocities as low 
as v out = 0.3 km s -1 , based on the spectral fit shown in 
Figure [To) we find that the total mass may increase by 
an order of magnitude or more, which was also shown by 


Down e s fc C abrit (2007); Offner et al. (2011); Dunham 
et al.| (|2014D lor low-velocity outflows. When assessing 


the impact that outflows have on the cloud environment, 
the momentum and energy are more appropriate diag¬ 
nostics than mass, as they will be less affected by the 
low-velocity channels (i.e., a factor of a few when extrap¬ 
olating based on Figure 10) so that the main conclusions 
that we propose below will hold qualitatively. 

Conversely, the high-velocity outflows contribute more 
significantly to the momentum and energy budgets, 
rather than mass. The velocity range that we probe here, 
with \v out \ < 31 km s _1 , is mostly representative of the 
SHV outflow component. However, EHV components 
have been shown Jo be driven by some low-mass proto¬ 
stars (e.g., Tafalla et al. 2004), and if present they can 
account for 2-4 times as much momentum and energy 
as the SHV components. In the following sections we 


compare the relative contributions of momentum and en¬ 
ergy from outflows, turbulence, and gravity of the cloud. 
Although the magnitude of the relative contribution by 
outflows may change if in fact there are EHV components 
associated with these sources, our discussion and conclu¬ 
sions (below) that outflows significantly impact the sur¬ 
rounding environment hold. Nonetheless, given that we 
do not probe the entire range of low- and high-velocity 
outflow gas, it should be kept in mind that we likely 
under-estimate the mass and subsequent m ass-depe ndent 
kinematic properties expressed in Sections [4 


4. DISCUSSION 

4.1. Characteristic outflow properties 

In order to characterize the outflows in this region, 
we present several common outflow properties in Table 
[4] We list the dependence of each property on assumed 
inclination in column 2, where £ is the angle between 
the outflow axis and the line of sight (i.e., £ = 90° cor¬ 
responds to an edge-on system). For each property, we 
present values corresponding to the lower and upper lim¬ 
its of the characteristic length (columns 3-4); we also 
calculate the values corresponding to the median charac¬ 
teristic length or dynamical time, both uncorrected for 
inclination angle (column 5) and assuming an average in¬ 
clination angle of (£) = 57.3° (column 6). The respective 
units are shown in column 7. 

The first outflow property is the characteristic out¬ 
flow length ( R c har,out), and this va lue propagates in 
subsequent calculations. In Section |3.2[ we discussed 
several coherent outflow features which have lengths of 
45" — 210" ( R c har,out = 0-1 pc to 0.4 pc, respectively), 
and a median of 90" (R c har,out = 0.19 pc, or 0.22 pc for 
£ = 57.3°). Based on the morphology of the emission 
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TABLE 4 

Outflow Properties 


Property 

Inclination 

Dependence 

Limits a 

Lower U pper 

Medians 13 

Uncorrected (£) = 57.3 U 

Units 

Characteristic length ( R c har,out ) 

V sin£ 

0.1 

0.4 

0.19 

0.22 

[pc] 

Characteristic velocity ( V c har,out ) 

1 / cos£ 



4.8 

8.9 

[km s *] 

Dynamical time ( t dyn ) 

cos £/ sin£ 

1.9 

8.9 

3.8 

2.4 

Sh 

>1 

o 

1—1 

2 ^ 

Mass loss rate (M) 

sin£/ cos£ 

1.5 

0.3 

0.7 

1.2 

[xlO -4 M© yr" 1 ] 

Outflow force ( F ou t ) 

sin£/ cos 2 £ 

7.2 

1.5 

3.6 

10.3 

[xlO -4 Mq km s _1 yr -1 ] 

Outflow luminosity ( L ou t ) 

sin£/ cos 3 £ 

42.2 

9.0 

21.1 

112.6 

[xl0- 2 L o ] 


a “Lower” and “upper” correspond to outflow property listed in first column obtained using the lower and upper limits of R c har 
(row 1), respectively, and uncorrected for inclination angle. 

b Correspond to the outflow properties listed in first column obtained using the median R c har value (row 1), with “uncorrected” 
indicating raw values without correcting for inclination angle, and (£) = 57.3° is the assumed average outflow inclination angle 
with respect to the line of sight. 


detected here, we suggest that the majority of the out¬ 
flows that we discussed originate from the YSOs near the 
dense cluster center in the map, and we assume that the 
previously mentioned values for R c har,out are the limits 
and median lengths representative of outflows found in 
this region. We use these values in further diagnostic cal¬ 
culations, although it has also been shown that outflow 
lengths var y by at least an order of magnitude in other re - 
gions (e.g., Nakamura et al.|2011b||Plunkett et al.|2013|), 
so more extreme values may be feasible. Nonetheless, 
our median value and range that we present in Table [4 
are consisten t with the literature, especially Naka mura 
et al. (2011a) who estimated the mean length of outflow 
lobes in Serpens South to be 0.1 — 0.2 pc, or 0.17 — 0.33 
pc at the assumed distance of 429 pc. 

Next, we calculated the characteristic velocity 
(V char,out) which is the total outflow momentum di¬ 
vided by outflow mass in the region, and we found that 
Vchar,out = Pout/M out = 4.8 km s _1 . Correcting for 
an inclination angle £ = 57.3° resulted in V c har,out = 
8.9 km s -1 . The resulting median dynamical time de¬ 
fined as t dyn = Rchar,out /v char , out is 3.8 xlO 4 yr (2.4 xio 4 
yr for £ = 57.3°), with a range from 1.9 — 8.9 x 10 4 yr for 
the Rchar,out range mentioned previously. The remain¬ 
der of diagnostics in Table 0] pertain to outflow mass, 
momentum, and energy, and we explain those diagnos¬ 
tics in the following sections. 


4.2. 


Outflow momentum and energy compared with 
turbulence 


In the model of Nakamura & Li (2007), outflows bal¬ 
ance infall motions and accretion, sustain the cloud in 
quasi-static equilibrium, and lead to a slow SFR. With¬ 
out outflows to inject additional turbulent motions in this 
model, the turbulence would dissipate in less than sev¬ 
eral free-fall times. Therefore, the following calculations 
aim to quantify momentum and energy injection by out¬ 
flows, specifically in the context of infall and turbulence 
in the same region. The calculations of mass loss rate, 
outflow force, and luminosity (see rows 4-6 of Table [4]) 
are inversely proportional to the dynamical time, which 
is a source of uncertainty of a factor of a few as explained 
in the previous section. We report the median values in 
the following text, but Table 0] also includes a range that 
corresponds to the lower and upper limits for t dyn (un¬ 
corrected for inclination angle), as well as median values, 


both uncorrected and corrected for an average inclination 
a ngle of (£) = 57.3°. 


Reipurth & Bally (2001) explained that the “unified” 
outflow model relates outflows with jets and/or winds, 
such that the same primary engine powers a YSO’s high- 
velocity jet/wind as well as lower-velocity lobes of an 
outflow traced by CO emission. Hence, the cumulative 
mass loss rate via outflows and the mass loss rate via 
winds can be expressed according to M owt V char,out = 

M-windVwind • From our CO observations, we found 
the cumulative mass loss rate via molecular outflows 
(Mout = Mout/tdyn) to be 1.2 x 1CT 4 M 0 yr -1 (assum- 
ing £ = 57.3°). It follows that the cumulative wind 
mass loss rate is M w i nd — 1-0 x 10 -5 M 0 yr -1 , as¬ 
suming V W i n d = 100 km s -1 (with _a like ly ran ge of 
V„rind = 10 0 — 200 km s _1 , c.f. Bally et al. 2007 Arce 
et alj2007 , and references therein). Further/it has been 


shown that mass loss rates are tied to accretion rates 

Based on the wind mass loss 


(Hartigan et al. 2000). 
rate, and assuming the fraction of stellar mass launched 
in a wind to be f w ^ 0.1 — 0.3 (|Shu et al.||1988l [Pelletier 
& Pudritz||1992| Lee et al.||2UU7a|b|), we expect the total 
mass accretion rate onto the stars 


this region to be ^ 0.3 — 1.0 x 10 4 M@ yr 1 (or lower 


powering outflows in 

© yr“ - 

by a factor of 2 if we assume the larger V w i nd = 200 
km_s -1 ). 

Kirk et al. (2013]) estimated the infall rate of material 


onto the cluster-forming clump to be about M in f a u 
2.2 x 10 -4 Mq yr -1 (assuming filament inclination an¬ 
gle that they find to be 12°, and scaling to a distance 
of 429 pc). We can relate this quantity with the to¬ 
tal accretion rate onto stars within the cluster using 
Mace = f sFEMjnfaii, with fs FE typically in the range 
of 0.1-0.3 (Lada & Lada p()03 ), since not all the material 
that is added to the clump through gravitational infall 
will be used to make protostars (and disks). Assuming 
f w = 0.1 and fsFE = 0.3, we then estimated the infall 
rate for this region to be ^ 3 x 10 -4 M 0 yr -1 , which 
is compar able to t he infall rate onto the cluster found 
by Kirk et al. (2013). Hence, the cumulative mass loss 
rate via molecular outflows that we measured and the 
observed infall rate appear to be consistent in Serpens 
South. 

In addition to assessing the outflowing mass over time, 
we calculated the total outflow momentum per unit 
time, known as the cumulative outflow force (F out = 
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Pout/tdyn ), to be 10.3 x 10 -4 M 0 km s _1 yr _1 (assuming 
£ = 57.3°). Further, the cumulative outflow lum inosity 
(L out = E nut /tdim) is 1.13 1/0. For comparison, Naka-|| 
mura et al. (2011a| reported a range L out = 0.8 — 1.7 L 0 
^Tbrrected for a distance of 429 pc), assuming an incli¬ 
nation angle of £ = 57.3° and a range of characteristic 
lengths R out = 0.1 — 0.2 pc. 

Comparing the turbulence energy dissipation rate and 
the outflow luminosity quantifies the contribution of out¬ 
flows and their impact on the cloud turbulence. In or¬ 
der to determine this, first we solved for turbulent en¬ 
ergy, which is defined as E tur b = \M c ig\ d . Through¬ 
out, we adopted a cloud mass of M c i = 220 M 0 and 
a radius of R c \ = 0. 17 pc, based on observ ations of 
HCO + ( J = 4 — 3) by Nakamura et al. (2011a) and ad¬ 
justing for a distance of 429 pc; the mass is consistent 
within ~ 20% of that estimated by Kirk et al. (2013) 
for the central cluster. Also, 3D and ID velocity dis- 
persion were related according to the following relations: 
(J 3 d = V3<j v and a v = A v tur b /(2 y/2 In 2). The veloc¬ 
ity dispersion A v tur b = 1-36 km s -1 (Fernandez-Lopez, 
2014, private communication) was determined from the 
FWHM of N 2 H+ which, traces mostly dense and qui¬ 
escent gas. The resulting turbulent energy within the 
mapped region was found to be 2.2 x 10 45 erg. 

The turbulence energy dissipation rate is given by 
Lfurb = E tur b/tdiss • The diss ipation time is defined as 
tdiss = tff (3.9 k/M r ms) (c.f. M ac Lowp9 99), where free- 
fall time is tff = y/3i r/(32 Gp) , and ft = A^/Aj ^ 1 as¬ 
suming that the turbulence driving length (A^) is equal 
to the Jean’s length (Aj) of the clump. We deter¬ 
mined the gas density of the cloud (p = 3 M c //47tR 3 j) 
to be p = 8 x 10 -19 g cm -3 . The Mach number is 
Mr ™,, a = as da/ 2.72m fir /(fcT) = 5.5, with T = 10.8 
K (Friesen et al. 2013). The resulting turbulent en¬ 
ergy dissipation rate is L tur b = 0.3 L 0 . The ratio 
tl = L ou t/L tur b ~ 0.6 — 3 (the range corresponds to 
L out uncorrected for inclination angle, and assuming 
£ = 57.3°) quantifies the potential for outflows to coun¬ 
teract turbulence dissipation in the region. A value of at 
least 1 indicates that outflows drive energy comparable 
to or greater than the turbulent energy in the region, as 
is the case here after correcting for inclination angle. 

Another way to assess the potential impact of outflows 
on the cloud turbulence is to compare the total outflow 
force and the turbulence momentum dissipation rate. In 
the outflow-regulated cluster formation model, the out¬ 
flow force should be comparable to or greater than the 
turbulence momentum dissipation rate in order for out¬ 
flow feedback to replenish internal supersonic turbulence 
and maintain the clump close to quasi-viria l equilibrium 
(|Li & Nakamura|2006||Nakamura & Li|2014|). Turbulence 
momentum dissipation rate of the internal turbulent mo¬ 
tion is d efined as dP tU rbldt = — 0.21 M r j/tdiss (equa¬ 
tion 4 of Nakamura & Li 2014| ). We found dP tU rb/dt — 
—2.8 x 10 -4 Mq km s -i yr -i , which is approximately 
25% of the median outflow force (assuming £ = 57.3°). 
Even accounting for a factor of a few uncertainty in 
the outflow force that we mention above, the ratio of 
dPturb/dt to F ou t likely remains less than 1. Assuming 
that outflows can transfer momentum and energy to a 
substantial fraction of the gas in the immediate vicin¬ 


ity of the dense cluster of young protostars, the ratio of 
outflow force to turbulence momentum dissipation along 
with the ratio vl (see above) suggest that outflows are 
an important source of turbulence in this region. 

4.3. Outflow energy compared with gravitational energy 

We compared the outflow energy with the gravitational 
binding energy of the clump to investigate whether the 
outflows could potentially unbind the cloud. Using the 
previously mentioned cloud mass and radius, the gravi¬ 
tational binding energy (W = —GM^/R c i) of the cen¬ 
tral region is —2.5 x 10 46 erg. The cloud mass is about 
50 times greater than the outflow mass, while the grav¬ 
itational binding energy is about 7 times greater than 
the outflow energy, assuming an inclination angle of 
£ = 57.3° (see Table [3]). In order to quantify whether 
outflows have sufficient energ y to overcome gravity and 
disperse the surrounding gas, Nakamura & Li| (|2014|) uti¬ 
lized the non-dimensional parameter r\ ou t = —2 E 0 ut/W. 
We found that r) ou t — 0.3 (assuming £ = 57.3°), which 
implies that outflow feedback is not sufficient to provide 
enough kinetic energy to unbind the entire cloud. 

4.4. Serpens South and NGC 1333: a comparison 

Serpens South is the second region ob served as part 
of a s tudy of clustered star formation; in Phm kett et al.| 
(2013) we mapped the prototypical protostellar cluster 
NGC 1333 utilizing CARMA and FCRAO observations 
of CO J = 1 — 0. Given the consistency among the 
observations, here we compare Serpens South and the 
more evolved NGC 1333, and we present an empirical 
scenario of outflow-cluster interaction at different evolu¬ 
tionary stages. The comparisons that we discuss below 
a re also summarized in Ta ble |U 


Gutermuth et al. (2008) suggested that the high per¬ 
centage of YSOs that are protostars is evidence of Ser¬ 
pens South’s youth. Here we present a tally of YSOs 
in Serpe ns South using the most recent classification s 
following Evans et al. (2007) and Maury et al. (2011), 
allowing a consistent comparison with 1NGG 13351 We 
count Class 0, Class I, and flat-SED sources as proto¬ 
stellar sources, and these along with Class II/III sources 
comprise YSOs. We count 45 protostellar sources and 
a total of 55 YSOs within the ~ 22 arcmin 2 (0.3 pc 2 ) 
region of Serpens South that we show outlined with a 
white rectangle in Figure [2j of which 32 protostars and a 
total of 37 YSOs reside within the ^ 8 arcmin 2 (0.1 pc 2 ) 
inner-most cluster region outlined by the red rectangle in 
Figure |2j We compare this to the 26 protostellar sources 
out of a total of 55 Y SOjMnt he 0.23 pc 2 region of NGC 
1333 mapped by Plunkett et al. (2013). 

The protostellar number densities (number of proto¬ 
stars per unit area mapped) are 130 and 110 pc -2 , re¬ 
spectively, for Serpens South and NGC 1333, showing 
that both of the regions we mapped are very densely 
populated with young members but about 20% more so 
in Serpens South. Further, in Serpens South the proto¬ 
stellar fraction (the ratio of protostellar sources to YSOs) 
is about 80%, whereas in NGC 1333 this percentage is 
50%. We also note that within the central region, out¬ 
lined by the red rectangle in Figure J2[b), the protostellar 
fraction of Serpens South rises to 90%; the central region 
corresponding to the same physical area in NGC 1333 
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TABLE 5 

Comparison of Serpens South and NGC 1333 


Property 

Serpens South a 

NGC 1333 a 

Units 

Protostellar number density 

130 

110 

pc -2 ] 

Protostellar fraction 

80 

50 

%] 

Outflow mass per area 

9.7 

26 

M© pc 2 ] 

Outflow emission prevalence 13 

90 

60 

%] 

r L (— L out /L turb ) 

0.6 - 3 C 

6 - 24 c 


'Hout (— 2 E/ on ^/fT) 

0.1 - 0.3 C 

0.2 - 0.7 C 



Pertains to area enclosed in white rectangle in Figurej^Jfor Serpens South, and region mapped 


by |Plunkett et al.||2013^ for NGC 1333. 

Outflow emission coverage (pixels projected onto the plane of the sky) with respect to a mapped 
area of 0.3 pc^. 

C The given lower- and upper-limits of ranges correspond to (1) raw values without correcting 
for inclination angle, and (2) assuming (£) = 57.3° as the average outflow inclination angle with 
respect to the line of sight, respectively. 


maintains a protostellar fraction of ^ 50%. The higher 
percentage of young sources in Serpens South suggests 
that the majority of sources in this region have only re¬ 
cently begun to drive outflows that impact the surround¬ 
ing cloud environm ent, compar ed with NGC 1333. 

In NGC 1333,|Plunkett et al.|(|2013|) identified outflows 
associated with 6 of 11 Class 0 sources and 2 of 11 Class 
I sources, but no outflows were associated with more 
evolved sources. This corresponds to 8 of 26 protostel¬ 
lar sources driving identifiable outflows, and at least four 
candidate outflows had unidentifiable driving sources. It 
is worth noting that NGC 1333 is clos er (235 pc) than 
Serpens South (429 pc, see Section 1.2), and its outflow 
features are more clearly distinguishable such that we 
were able to associate outflows with their driving sources. 
This is not the case for Serpens South, which is farther 
and has a much more complex web of interspersed out¬ 
flows, as seen in Figure [3] and elsewhere, so our compar¬ 
isons her e are in terms o f cum ulative outflow activity. 
In total, Plunkett et al.| ([2013) measured 6 M 0 of gas 
associated with outflow lobes or outflow candidates in 
NGC 1333 within a region of ~ 0.23 pc 2 , correspond¬ 
ing to ^ 26 M 0 pc -2 . In Serpens South, we measure 
~ 2.8 M 0 of CO gas within a region of ^ 0.34 pc 2 , or 
~ 10 Mq pc -2 . Granted, a comparison of mass entrained 
by outflows relies on the assumption that the initial con¬ 
ditions of the cluster before star formation began are also 
comparable. 

Spatially, outflow emission is very prevalent in both re¬ 
gions. To quantify the fraction of the mapped area with 
outflow emission, we measured how many pixels (pro¬ 
jected onto the plane of the sky) have 12 CO signal-to- 
noise ratio greater than three in at least three velocity 
channels (“outflow emission prevalence” in Table [5]). In 
the Serpens South map, 80% of the mapped region has 
emission in at least three velocity channels. For com¬ 
parison, we scaled em ission in the NGC 1333 map from 
Plunkett et al. (2013) to a distance of 429 pc and the 
same rms noise threshold as the Serpens South data. In 
a region corresponding to the same physical area of the 
Serpens South map we measured emission (again in at 
least three velocity channels) spanning about 60% of the 
mapped area on the plane of the sky. The most simple 
explanation for this is that the number density of pro¬ 
tostellar sources, assumed to drive outflows, is higher in 
the central-most region of Serpens South than in NGC 


1333. A more detailed study of individual outflows in 
Serpens South is needed to develop this argument. Con¬ 
sidering that NGC 1333 is a more evolved protostellar 
region, we expect some of the outflowing material to 
have escaped that region, whereas in Serpens South the 
material has had less time to travel beyond the central 
cluster that we mapped. This interpretation of outflow 
emission prevalence compared between the two regions 
assumes that outflows are randomly oriented in both 
regions. Alternatively, if large-scale filaments and gas 
kinematics cause preferential orientations of the outflows 
and/or their driving sources, then the spatial prevalence 
projected on the plane of the sky would depend on the 
inclination angle of the majority of outflows. Since Ser¬ 
pens South is known to have a concen trati on of YSOs 
along a dense filamentary structure (Gutermuth et al. 


2008), then this structure potentially affects the spatial 
distribution of outflow emission. 

Two previously mentioned ratios, vl and r] outl quantify 
the relative contributions of energies from outflows and 
turbulence, and from outflows and gravity (respectively) 
in these regions. For the slightly more evolved region 
NGC 1333, these parameters are Cl ~ 6 — 24 and r\ out = 
0.2 — 0.7 (for £ = 0 — 57.3°), whereas in Serpens South 
we found these values to be vl ~ 0.6 — 3 and rj out = 
0.1 — 0.3 (for £ = 0 — 57.3°), respectively, which are 
about 10-50% the values for NGC 1333. The outflows in 
both regions contribute sufficiently to drive turbulence, 
although to a greater extent in NGC 1333, and in NGC 
1333 the outflows may nearly counter the gravitational 
energy of the clump. It follows that as a cluster evolves 
from the very young stage of Serpens South to the slightly 
more evolved stage of NGC 1333, the outflows sustain 
turbulence in a similar proportion, but destruction via 
outflows is more imminent as a region evolves. 

Finally, we explore an evolutionary scenario in which 
Serpens South may come to resemble NGC 1333. First, 
we calculated the relative protostellar fraction within 
Serpens South as it pr ogresses f rom its current, young 
age (~ 2 ± 1 x 10 5 yr, F riesen et al.||2 013) to the slightly 


more evolved state of NGC 1333 over its expected age of 
1 — 2 Myr (|Lada et al.||1996). We assumed lifetimes 
of 0.10, 0.34, and 0.35 IVlyr spent as Class 0, Class I , 


and Flat SED classes, respectively (Evans et al. 2009| 
and a uniform distribution of ages within each of these 
classes. In the simplest (unlikely) case, if no additional 
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protostars form, then the protostellar fraction will reach 
that of NGC 1333 50%) in 0.4 Myr. In order to ac¬ 

count for the ongoing formation of additional, young pro- 
tostars, then incorporatin g also the star formation rate 
(SFR) of 90 M© Myr -1 (Gutermut h et al. 20081), and 
an average YSO mass of 0.5 AZ© (c.f. Evans et al. 2009), 
then we add one more protostar every 0.006 Myr, and 
in that case the protostellar fraction reaches 50% in 1.3 
Myr. However, the SFR in the more evolved region NGC 
1333 is about half of that in Serpens South, suggesting 
that the SFR likely decreases as a cluster evolves. If in¬ 
stead we assume the SFR over t he e ntire time to be 45 
M© Myr -1 (as for NGC 1333, Lada et al.|l996), then the 
protostellar fraction reaches 50% in 0.9 Myr. Additional 
work should be done to determine how the SFR changes 
over time for a cluster, but even based on these simple 
calculations using representative values for the SFR, the 
protostellar fraction of Serpens South will likely decrease 
to be comparable to that of NGC 1333 in ^ 1 Myr. 

Also, as Serpens South evolves to the age of NGC 
1333 during the next ~ 1 Myr, we expect mass to es¬ 
cape from the cloud, possibly also related to the re¬ 
spective power-law slopes of the mass spectra in Figure 
[To] We calculated the escape velocity of Serpens South 
to be v esc = 3.4 km s _1 , which is higher than that of 
NGC 1333 where v esc = 2.2 km s -1 (jArce et al. 2010), 
but nonetheless most of the outflow mass we calculate 
in both regions pertains to gas with outflow velocities 
greater than the escape velocity. Assuming ages of Ser¬ 
pens South and NGC 1333 as above, and assuming a 
constant outflow rate, then the current mass loss rate 
of 70 — 120 M 0 Myr -1 implies that ~ 100 M© may be 
lost before Serpens South reaches the evolutionary stage 
of NGC 1333. At that time, with a cloud mass that is 
about half of what we assumed here for the calculations 
of gravitational energy of the clump, the value of rj out 
in Serpens South will approximate that of NGC 1333. If 
cumulative outflow energy decreases with time, then Ser¬ 
pens South would come to resemble NGC 1333 sooner. 
Although this is a simple order of magnitude calculation 
with admittedly simple assumptions, we consider it valu¬ 
able to place Serpens South and NGC 1333 within a rel¬ 
ative evolutionary scenario. More thorough treatment of 
the assumpti ons we made here, including accretion onto 
the filament (|Kirk et al.||2013) and varying (and overall 
decreasing) mass ejection rate over time, as well as fur¬ 
ther detailed analysis of additional clusters spanning a 
range of evolutionary stages will constrain this scenario. 


5. SUMMARY 

Here we assess the activity of outflows in the young pro¬ 
tostellar region Serpens South, specifically focusing on 
the central 0.3 pc 2 where we find an intricate web of out¬ 
flow activity. We present observations of CO J = 1—0 for 
which interferometer (CARMA) and single dish (IRAM 
30 m telescope) observations were combined to probe a 
range of spatial scales from approximately 0.01 pc (or 
2000 AU) to 0.8 pc (at a distance of 429 pc). These 
provide the highest-resolution observations of molecular 
outflows in this region to date, while subsequently retain¬ 
ing extended cloud-scale emission. 

Utilizing the CO J = 1 — 0 data from CARMA and 
the IRAM 30 m telescope, in addition to CO J = 3 — 2 
data from APEX and CSO, we determined an opacity 


correction factor for the main beam temperature of the 
12 CO observations as well as solved for varying excitation 
temperature throughout the region; these methods are 
described in Appendix [A] and [B] These observations and 
methods were necessary to obtain accurate estimates of 
outflow mass, momentum and energy. The results can 
be summarized as the following: 

1. Within the region we mapped, we measured the 
total outflow mass to be M out = 2.8 M©, outflow 
momentum to be P out = 13.6 M© km s _1 , and en¬ 
ergy to be E out ^ 1.0 x 10 45 erg (uncorrected for 
assumed inclination angles of outflows). Several 
methods using different assumed excitation tem¬ 
peratures resulted in masses, momenta and ener¬ 
gies that varied by factors of a few. 

2. The outflow mass spectra were fit with power laws 
of 3.5 zb 0.1 and 3.3 ± 0.1 for the blueredshifted 
and redshifted emission, respectively, for velocity 
channels with 3 < \v out \ < 15 km s _1 . We found 
steeper slopes for the outflows in the more evolved 
region NGC 1333 than for Serpens South, and we 
found no break for either region, having properly 
accounted for optical depth effects and recovered 
the “hidden” mass at low outflow velocities. Likely 
a shallower (steeper) slope is a sign of the earlier 
(later) evolutionary state of a cluster. 

3. We quantified the relative contributions of outflows 
compared with turbulence and gravity in this re¬ 
gion, comparing the outflow energy input rate to 
the turbulence energy dissipation rate and the out¬ 
flow energy to the gravitational energy, as well as 
outflow momentum force to the turbulence momen¬ 
tum dissipation rate. We found that outflows con¬ 
tribute enough momentum and energy to maintain 
turbulence, but outflow feedback plays a minor role 
compared with the gravitational potential energy 
of the clump. These mass, momentum, and energy 
diagnostics are consistent with the quasi-static star 
formation model. Additionally, based on the cumu¬ 
lative mass loss rate, we estimated a mass accretion 
rate that is comparable to the lower limit infall rate 
onto the filament. 

4. In Serpens South, the bulk of the outflow emis¬ 
sion appears to be emanating from sources located 
within a central 0.1 pc 2 region of high protostellar 
density. Our 2.7 mm continuum observations with 
CARMA show seven distinct continuum emission 
sources that are considered protostellar candidates, 
some of which correspond to the ^ 45 Spitzer pro¬ 
tostars in this region. The high resolution of our 
observations reveals the extended morphologies of 
CARMA-3a/b and CARMA-6/7, and we measure 
projected distance separations for each proposed 
binary (or multiple) system of about 3000 and 6500 
AU, respectively. 

5. In a census of YSOs in Serpens South compared 
with NGC 1333, the central hub of the Serpens 
South cluster contains a relatively high protostellar 
fraction and protostellar number density, and we 
compare outflow energetics and their influence on 
the respective clouds. 
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Considering these outflow diagnostics and the YSO 
census in Serpens South, we propose an empirical evo¬ 
lutionary scenario for which Serpens South reveals an 
early stage of star formation, and this complex region 
will likely evolve to resemble a protostellar region like 
NGC 1333 within approximately 1 Myr. Even higher 
resolution continuum emission observations are needed to 
reveal several sources as likely multiple-member systems, 
and higher sensitivity observations will reveal weaker 
sources, refining our view of the young sources that 
drive outflows in this region. Similarly, higher resolu¬ 
tion molecular outflow observations, of which ALMA is 
capable (A.L. Plunkett et al. in prep), will allow us to as¬ 
sociate specific outflows with their driving sources, and 
better characterize individual outflows as was done for 


NGC 1333 by Plunkett et al. (2013). 
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APPENDIX 


CORRECTING FOR OPACITY 


In order to accurately measure masses and energetics of outflows, Dunham et al. (2014| argue that careful attention 
must be paid to account for the effects of line opacity and excitati onTemperature variations (see Appendix IbI), among 
other factors. We follow the method presented by Dunham et al. (2014, and references therein) to correct for opacity, 
given that low-J transitions of 12 CO are typically optically thick at cloud systemic velocities, and this likely needs to 
be accounted for in the line-wings as well. Here, opacity is determined for each velocity channel, independently for the 
J = 1 — 0 and J = 3 — 2 energy level transitions. For the J = 1 — 0 transition, we used the combined CARMA+IRAM 
maps, and the 13 CO data cube was regridded to have the same pixel size and channel width as the 12 CO data cube. 
Units of the CARMA+IRAM maps were converted from Jy beam -1 to Kelvin by dividing by 3.78 and 4.13 for 12 CO 
and 13 CO, respectively. For the J = 3 — 2 transition, we used data from APEX and CSO for 12 CO and 13 CO, 
respectively, after convolving the APEX data cubes to have the same beamsize (35") as the CSO observations, and we 
smo oth ed th e CSO spectra to have the same velocity resolution (0.2 km s -1 ) as the APEX observations (see Sections 


2.3 


and 


2.4) 


Figure 11 shows the channel-by-channel intensity of 12 CO and 13 CO. For the J = 1 — 0 maps, we averaged within 


a box 60 x 60" centered at R.A.= 18 /l 30 m 03 s , decl.= —02°02 / 57", where a high prevalence of outflow emission is 
observed. In the case of the J = 3 — 2 maps, we averaged the region of each map corresponding to the two CSO 
pointings nearest to the center (see Figure HI). Independently for J = 3 — 2 and J = 1 — 0, we fit a parabolic function 
to the ratio Ri 2 /i 3 (^) = /(Ii 3 (vi)), where (Ii 2 {vi)} and (/i 3 (t^)) are averaged intensities for each velocity 

channel Vi. The parabolic fits are described with the following equations: 


J = 1-0: R 12 /is(vi) = (1-75 ± 0.03) (v* - 8.14 ± 0 . 02) 2 + (-3.2 ± 3.2) (Al) 

J = 3 - 2 : R 12 /i3(vi) = (1.19 ± 0.01)(v* - 8.78 ± 0.01) 2 + (1.6 ± 1.6). (A2) 


We fit the function to the velocity ranges (selected by eye to reasonably fit a parabolic function) indicated with red 
points in Figure [TTJ and we truncated each function with a maximum value of 62 at the high-velocity edges of the 
spectra. Also, we do not constrain the fits to reach a minimum at the cloud velocity, as can be seen from the equations 
above which both have minima at velocities greater than 8 km s -1 . This asymmetrical opacity correction may be an 
effect of averaging over a region that is comprised of several outflows with distinct orientations and inhomogeneous 
initial cloud gas distribution, and therefore the averaged outflow emission is inherently asymmetric. Using these 
functions, we numerically determined the velocity-dependent optical depth t\ 2 for both transitions, following the 
relation: 
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Fig. 11. — Top row shows 12 CO (solid green line) and 13 CO (dashed magenta line) spectra for the J = 1 — 0 (left) and J = 3 — 2 
(right) transitions, averaging over the regions specified in the text (see Appendix |A|. Bottom row shows the ratio of these spectra for both 
transitions, R 12/13 — -^ 12 /^ 13 - A parabola, shown with a red dotted line, was fiFTo the points also shown in red for each transition, and 
the equations for these parabolas are given in Equations and ( |A2| , truncated at a value of 62. Channels with 13 CO emission below 
3cr and/or i?i 2 /i 3 > 62 are excluded from the plot. 


R 


12/13 


[ 12 CO] /l-exp[—r 12 ] 
[ 13 CO] V r 12 


(A3) 


assuming that the 13 CO is optically thin. We then correct the observed main-beam temperature, 12 , for the 12 CO 


velocity-dependent opacity, resulting in the corrected main-beam temperature, T r 


mb, 12’ 


given by: 


T '* _ rp 

mb, 12 - -‘-mb, 12 


Tl2 


1 - exp(—ri 2 )' 


(A4) 


In Figure 12 we show the averaged 12 CO spectrum versus velocity before and after making this opacity correction, as 


well as theTactor that relates observed main-beam temperature and corrected main-beam temperature. The correction 
factor shown in Figure [l2|b) reaches ' 
velocities with \v LS r - v doud \ >6kms 


5 for velocity channels with |vlsr — v c i ou d\ 
-1 


3 km s 1 and reduces to 1 for 


EXCITATION TEMPERATURE 


In addition to opacity c orrectio n , excita tion temperature (Tex) variations have an effect on measuring outflow 
characteristics. Specifically, Dunham et al. (2014) show that outflow mass is uncertain by a factor of 0.7-3 depending 
on assumed T e x in the range ot 10-200 K. Excitation tempera tur e is g enerally assumed to be within the range of about 
10 - 50 K and perhaps higher (see discussion in Dunham et al. 2014) in protostellar regions such as Serpens South. 

Given the range of plausible values, rather than adopt “single value for the entire region, we opted to obtain the value 
on a pixel-by-pixel + channel-by-channel basis using the 12 CO J = 1 — 0 and J = 3 — 2 data. We used observations 
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Fig. 12.— (a) Average 12 CO J = 1 — 0 spectrum within the central ~ 16 arcmin 2 region of the map. The green (solid) line shows the 
12 CO emission without any co rrect ion for opacity, and the magenta (dashed) line shows the spectrum after correcting for optical depth of 
the line, according to Equation | A4| Vertical blue and red dotted lines mark the velocities v c i OU( i =b 3 km s — 1 , and we only consider velocity 
channels outside of this ra nge to be mostly unaffected by cloud emission, (b) Factor used to correct 12 CO main beam temperature for 
optical depth (see equation |A4|. 



Fig. 13.— Ratio of intensities I 10 /I 32 expected as a function of excitation temperature. See Appendix [B] where we describe how we 
account for varying Tex in our map. The vertical dashed green line marks the imposed minimum at Tex = 1TTK. Several sample error bars 
show uncertainty in Tex, based on the rms of the 12 CO J = 3 — 2 and J = 1 — 0 single-dish maps. 


from IRAM ( 12 CO J = 1 — 0) and APEX ( 12 CO J = 3 — 2 ) to solve for H 2 column density, 7V# 2 , as a function of 
excitation temperature according to equations [ 2 ] and [3j We did so independently using single dish observations of 
these two transitions, and then set the two equal to solve for T e x in each pixel and velocity channel. In other words, 
excitation temperature can be determined from the ratio of J = 1 — 0 to J = 3 — 2 intensities, I 10 /I 321 as we show in 
Figure [13] and according to the following equation: 


(Ej+i, 


10 


E 


J+1,32 ) 


In 


^10 (J10 + 1 ) I32 
_^32 (J 32 + 1 ) ho 


-1 


(Bl) 


where we used the opacity-corrected main-beam temperature (T^ b 12 , see equation 
ho and I32 

erg. Similarly, for the J = 3 — 2 transition we use the following values: J 32 = 2 , Ej+ 1^2 = 4.58 x 10 -15 erg, and 


A4) for each 12 CO transition as 

respectively. We used the following values for the J = 1 — 0 transition: Ji 0 = 0, Ej+ i ? i 0 = 7.64 x 10 -16 
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Ej+ lt32 /k = 33.19 K. 

Solving for T ex in pixels where iio and I 32 are above 4a, we measured excitation temperatures to be in the range 10 
K < T e x < 100 K. Where we could not reliably determine T ex for a given pixel, if we could do so for neighboring pixels, 
we took the averaged T ex in the surrounding box of 5 x 5 pixels and velocity channels within v out = ±0.5 km s -1 . 
Finally, we set a lower limit of T ex = 10 K for any pixels that did not meet the 4cr criterion, as well as pixels beyond 
the boundary of the ~ 5' x 5' APEX map (see Figure [I]). 

We then regridded the excitation temperature data to the grid of the CARMA±FCRAO map so that each pixel in the 
combined interferometer and single dish map, used for determining mass, had a corresponding excitation temperature. 
The map of excitation temperature is shown in Figure [15] The pixel-by-pixel ± channel-by-channel method described 


here, which we call method 1, results in an ave rage excitation te mperature of T ex (x, y, v) = 17 ± 8 K. For comparison 
with another protostellar region, in NGC 1333 (|F4unk ett et al.| 2013) we performed the same velocity-dependent pixel- 
by-pixel approach and found T ex (x,y,v) = 20 ± 4 K. In Serpens South, we also calculated T ex (v) channel-by-channel 
(method 2), where T ex (v) is the excitation temperature obtained for each channel by spatially averaging over the 
region (again, using the opacity-corrected main-beam temperature) and adopted for all pixels in that channel. By 
method 2, we found T ex (v) = 15 ± 4 K. Finally, we present three methods (3-5) that assumed constant values of 
Tex = 10, 30, 50 K to investigate this effect on mass. The distributions of excitation temperature values in the 
pixel-by-pixel ± channel-by-channel (method 1) and channel-by-channel (method 2) approaches are shown in Figure 
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Fig. 14. — (a) Normalized distribution of excitation temperatures determined pixel-by-pixel and channel-by-channel (method 1, magenta 
solid line) and spatially averaging Tex channel-by-channel (method 2, red dotted line). The vertical dashed green line marks the imposed 
minimum at Tex = 10 K. (b) Excitation temperature averaged for each velocity channel (method 2). The vertical dotted green line 
corresponds to the central cloud velocity of 8 km s -1 . 


We calculated the error in T ex by propagating the rms errors of /10 and I 32 in Equation [Bl 
(for method 1) incurred for several representative excitation temperatures in Figure 


13 


_ and we show the errors 
e error increases for higher 

Tex pixels; for example the errors in T ex are about ±1, ±10, and ±30 K for T ex =ll), 30, and 50 K, respectively. 

Mass, momentum and energy that we mention throughout the text (unless explicitly noted otherwise) are those 
determined using T ex (x,y,v) of method 1, taking into account the spatial- and velocity-dependence of excitation 
temperature, and in Table [3] we present these values for all T ex methods. Our five methods of T ex result in a factor 
of three difference for each outflow property (mass, momentum, and energy), with a constant T e x=50 K resulting in 
the highest values and constant T e x=10 K resulting in the lowest values (see Table |3|. Considering that the velocity- 
dependent pixel-by-pixel method shows that T ex ranges from 10 to 100 K within the region, then assuming a constant 
Tex will likely under-estimate (over-estimate) the mass in regio ns where the act ual T ex is higher (lower) than the 
assumed constant T ex . This is consistent with the work of Dunham et al. (2014), who showed that outflow masses 
determined using 12 CO J = 2 — 1 and 12 COJ = 3 — 2 varied by factors of 0.7 — 3 times depending on assumed excitation 
temperature in the range T e x=10-200 K, relative to the mass determined using their assumed T e x=50 K. 

Ideally, we would use several other energy level transitions in addition to the two that we used here to more accurately 
constrain the excitation temperatures as a function of position and/or velocity. Observations of high-J CO lines in 
other regions have revealed warmer excitation temperatures in the range 50 < T < 200 K associated with higher 


velocity gas (van Kempen et al. 2009a|b Yildiz et al.||2013| ). Likely the warmer gas is due to recent s hocks, while 


colder (and slower) gas is entrained after being shocked in the past. Based on higher-J transitions of CO, Yildiz et al. 
(2013|) measured higher excitation temperatures ( T ex > 50 K) than the mean values measured in our study, and they 
showed T e x to increase with outflow velocity. Hence, the high-J CO emission is important for probing the warmer 
and/or denser gas components, which are associated with broader line wings. 
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Our lower- J CO data correspond to the cooler, less dense gas components with narrower line wings; we show that 
excitat ion temperatures peak at an outflow velocity of about 5-10 km s -1 , and decrease at higher velocities. Since 


Yildiz et al. (2013) claim that the broad and narrow components contribute roughly equally, and we are only probing 
the cooler, narrow component, then only including the J = 1 — 0 and J = 3 — 2 transitions in our calculations likely 
underestimates the T e x measurement and subsequent mass, momentum, and energy determinations. It follows that 
the diagnostics th at we co nsider to assess the contribution of outflows to the energetics of the cloud, including tl and 
Vout (see Sections 4.2|4.3), should also be considered lower limits. Hence, the contribution of outflows may be more 
significant than we report here, reinforcing the main finding that outflows are important to the kinematics of this 
young cluster. 
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Fig. 15. — Maps of Tex for each pixel position, averaged at each pixel position within the channel intervals shown for each panel (these 
are the same intervals used in Figure^). Color scale bar is shown in the upper right, and chosen to truncate at T e x = 70 K to show color 
contrast in the maps. We measured excitation temperatures to be in the range 10 K < Tex < 100 K. 
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